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PLASMA PRODUCTION BY TRAVELING RESONANT PERTURBATIONS 


William I. Linlor 
Hughes Research Laboratories, Malibu, California 
(Received June 16, 1961) 


A major problem in fusion research is the pro- 
duction of plasma within a confining magnetic con- 
figuration. Attempts to ionize a gas initially pres- 
ent and to heat the ensuing plasma (as in the stel- 
larator machine) have been beset by difficulties,’ 
some of which have been explained theoretically.?»* 
Alternatively, attempts to inject energetic ions 
and electrons into mirror machines have been 
limited to low plasma densities when adiabatic 
conditions prevail, in agreement with theoretical 
analysis.*~® 

Nonadiabatic injection methods have been in- 
vestigated, such as change in charge-state” and 
entropy trapping.® Resonance trapping, for which 
the magnetic moment of a charged particle is in- 
creased by periodic interaction with magnetic field 
perturbations, has recently been proposed and 
tested, °~# 

In this note, we propose a new method of plasma 
production which involves traveling resonant per- 
turbations in an axially symmetric system. Axial 
symmetry is important, as otherwise loss of plas- 
ma would occur by drift motion across lines of 
magnetic flux. In order to avoid particle loss along 
paths that are mirror images of those which pro- 
duced trapping initially, a nonreciprocal effect 
must be present. This is provided by the use of 
traveling perturbations, which can be achieved by 
well-known techniques of spatial current distri- 
bution and polyphase excitation. 

Consider a modified mirror machine having a 
spatially periodic magnetic field on the axis of 
the form 


BoB, +B, cosan( f-2 “ved )s (1) 


where B, is the main magnetic field and varies 
slowly with time, the second term is a perturbing 
magnetic field having coordinate and time depen- 
dence as indicated, z is axial distance, \ isa 
wavelength and is a function of z, vy¢ is radio 
frequency, and? is time. The cyclotron frequency 
v, of a charged particle is given by 


os Be /2mm, (2) 


where e is the electric charge in emu, m is the 
mass in grams, and B is in gauss. The ratio 
vrf/Vvc is assumed to be essentially constant and 
much less than unity, typically of the order of 
0.01; the same is true for the ratio Bp/B ; 

The magnetic field specified by (1) is shown 
schematically in Fig. 1 at an arbitrary instant 
of time; the perturbation pattern moves from 
left to right at a velocity vyg measured in the 
laboratory frame. The wavelength \ is related 
to the z component vz of the velocity of the parti- 
cle in the laboratory frame by the relation 


Atv + Vg =e,- (3) 
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FIG. 1. Spatial variation of magnetic field intensity 
along axis. 
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It should be noted that as vz decreases, the value 
of decreases correspondingly. The decrease in 
vz is produced by conversion of velocity parallel 
to magnetic lines into perpendicular velocity. 
Transfer of kinetic energy between the particle 
and the traveling magnetic perturbation averages 
to zero because the particle is not in synchronism 
with the wave. During a given interaction cycle, 
only the “perpendicular” energy of the particle is 
slightly altered, but this does not affect the z com- 
ponent of the velocity of the particle. 

Because of the large difference in ion and elec- 
tron masses, simultaneous injection and trapping 
of a variety of positive and negative charges can 
be achieved in the same structure, provided that 
the appropriate resonance conditions for each 
species are satisfied. 

To describe the trapping mechanism for a given 
particle, let us consider an off-axis particle having 
an initial magnetic moment y, which enters at the 
left side of the magnetic field depicted in Fig. 1. 
In traversing the perturbations, the particle en- 
counters radial components of magnetic flux. As 
in the system of Sinel’nikov et al.,° azimuthal im- 
pulses are produced by the interaction of the longi- 
tudinal component of velocity with the successive 
regions of radial flux. These impulses, being in 
resonance with the cyclotron azimuthal motion by 
virtue of (3), cause the magnetic moment of the 
particle to increase to the value , while it passes 
from left to right, provided that the particle enters 
at the appropriate phase with respect to the rf 
perturbation. Upon reaching the right side of the 
system, the particle is reflected by the mirror 
field By- 

For particle motion in the right-to-left direc- 
tion, the resonance condition is far from satis- 
fied in this system, regardless of the phase of 
the azimuthal motion. Hence, the conversion of 
magnetic moment from an initial value to a larger 
final value is a “unidirectional” process for the 
system described. 

Subsequent traversals of the system by the now- 
trapped particle will cause the magnetic moment 
to fluctuate, >’? and eventually to decrease to a 
sufficiently low value that the particle is no longer 
contained by the mirrors, in accordance with 
Liouville’s theorem, unless time-dependent 
changes occur in the system. Such a change can 
be obtained by the adiabatic increase of B, with 
time, which permits filling phase space with in- 
jected and trapped particles. 

Evidently, the magnitude of the fluctuations of 
the magnetic moment of the trapped particles 
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depends sensitively on the resonance conditions; 
it is highly desirable that at least a hundred cy- 
clotron periods be required for conversion of the 
initial longitudinal velocity into perpendicular. 
These can be achieved by multiple passes within 
the same structure, employing superposed sharp- 
ly resonant systems. 

Interactions of the trapped particles with one 
another and with neutral particles, such as dis- 
sociation of molecular ions by collisions, scat- 
tering, etc., can produce sufficiently large changes 
in the resonance conditions that periodic decreases 
of Bz can be permitted without thereby producing 
loss of trapped particles. In such a system a 
practical Liouville limit in density is expected to 
be much greater than 10*° particles per cm’, which 
is a density of interest in fusion power production. 

The question may arise whether the traveling 
magnetic field would penetrate into the plasma 
when the latter has an appreciable conductivity. 
Skin depth has meaning only for conductors having 
fixed boundaries. For plasmas which undergo 
adiabatic magnetic field compression, there is no 
shielding of the interior because the guiding cen- 
ters of all of the particles loop constant flux @ 
(Busch’s theorem). Inasmuch as d¢/dt=0 for any 
selected element of plasma, no induced currents 
flow. In effect, during an adiabatic increase in 
the current in the external coils, the magnetic 
field within the plasma is compressed as if the 
plasma had infinite skin depth. 

A detailed analysis and description of the pro- 
posed system will be given in a later publication. 
A remark concerning phase stability of particles 
near resonance seems appropriate here. In Fig. 
2 we see that a reference particle having a phase 
angle a and exact resonance will have successive 
phase points labeled y. Other particles starting 
at y and having velocities larger or smaller than 
that of the reference particle will move to phase 
points / or s, respectively. However, a particle 
at 1 experiences a greater retardation than the 
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particle at ry; similarly a particle at s experiences 
a lesser retardation than the particle at r. Thus 
on successive cycles the points / and s move to 

l’ and s’, and thus closer to r, giving phase sta- 
bility. Obviously the argument is valid for -7/2 
<a<7/2. 

It should also be remarked that the proposed 
system described thus far has the objective of 
producing a “starting plasma” which subse- 
quently can serve as a dissociation agent for en- 
ergetic injected molecular ions. Once the start- 
ing plasma has been produced regardless of the 
complexity and cost, subsequent continuous feed- 
ing of fusible material seems to be relatively sim- 
ple. 
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GENERATION OF OPTICAL HARMONICS* 


P. A. Franken, A. E. Hill, C. W. Peters, and G. Weinreich 
The Harrison M. Randall Laboratory of Physics, The University of Michigan, Ann Arbor, Michigan 
(Received July 21, 1961) 


The development of pulsed ruby optical masers':” 
has made possible the production of monochro- 
matic (6943 A) light beams which, when focussed, 
exhibit electric fields of the order of 10° volts/cm. 
The possibility of exploiting this extraordinary 
intensity for the production of optical harmonics 
from suitable nonlinear materials is most appeal- 
ing. In this Letter we present a brief discussion 
of the requisite analysis and a description of ex- 
periments in which we have observed the second 
harmonic (at ~ 3472 A) produced upon projection 
of an intense beam of 6943A light through crystal- 
line quartz. 

A suitable material for the production of optical 
harmonics must have a nonlinear dielectric coef- 
ficient and be transparent to both the fundamental 
optical frequency and the desired overtones. Since 
all dielectrics are nonlinear in high enough fields, 
this suggests the feasibility of utilizing materials 
such as quartz and glass. The dependence of polar- 
ization of a dielectric upon electric field E may be 
expressed schematically by 


2 
paye(+E+ Eis), (1) 
where E,, E,... are of the order of magnitude of 
atomic electric fields (~10°® esu). If E is sinusoi- 
dal in time, the presence in Eq. (1) of terms of 
quadratic or higher degree will result in P con- 
taining harmonics of the fundamental frequency. 
Direct-current polarizations should accompany 
the even harmonics. 

Let p be that part of P which is quadratic in E; 
that is, p is a linear function of the components 
of the symmetric tensor EE. The eighteen coef- 
ficients which occur in this function are subject 
to restrictions due to the point symmetry of the 
medium. These restrictions are, in fact, identi- 
cal with those governing the piezoelectric coef- 
ficients. In particular, p necessarily vanishes 
in a material such as glass which is isotropic or 
contains a center of inversion. For crystalline 
quartz, however, there are two independent co- 
efficients a and 8 in terms of which 


‘ 2.672 
e. a By )+6E E 
py = “AEE - aaB B, 
p_=0 (2) 
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Table I. The square of the total p perpendicular to 
the direction of propagation of light through crystal- 
line quartz. 





Direction of incident The square of the total p 
beam perpendicular to direction 
of propagation 





= 2 $.. 
x (E. 0) p  * 0 
, = 2 2-2 4 
y G, 0) p, a a ES 

m 2 Sache 2 2)2 
z (E, 0) b, 0, a (EY +s, ) 





(z is the threefold, or optic, axis; x a twofold 
axis). If a light beam traverses quartz in one of 
the three principal directions, Eqs. (2) predict 
the results summarized in Table I. The second- 
harmonic light should be absent in the first case, 
dependent upon incident polarization in the second 
case, and independent of this polarization in the 
third. 

If an intense beam of monochromatic light is 
focussed into a region of volume V, there should 
occur an intensity J of second harmonic given 
(in Gaussian units) by 


I~ (w*/c*) (pu(V/v), (3) 


where w is the angular frequency of the second 
harmonic, c the velocity of light, and v an ef- 
fective “volume of coherence”; that is, the size 
of a region within the sample in which there is 
phase coherence of the p excitation. (This volume 
may in practice be much smaller than V.) An 
estimate of v is governed by several considera- 
tions. For example, it is probably of no greater 
extent in the propagation direction than ~ [n, 

x(n, -n,)"),, where n, and n, are the indices 

of refraction for the fundamental and second har- 
monic frequencies, respectively, and A, is the 
wavelength of the second harmonic. The lateral 
extent of this volume is determined in large part 
by the coherence characteristics of the optical 
maser. The situation for a maser of the gas 
discharge® type is clearly more favorable in 
this respect than that for the ruby device.’»? For 
a coherence volume of 10“ cm’, which we think 
may be realistic in our case, Eq. (1) indicates 
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FIG, 1. A direct reproduction of the first plate in which there was an indication of second harmonic. The 
wavelength scale is in units of 100 A. The arrow at 3472 A indicates the small but dense image produced by the 
second harmonic. The image of the primary beam at 6943 A is very large due to halation. 


that second harmonic intensities as high as a 
fraction of a percent of the fundamental could be 
achieved. 

In the experiments we have used a commerci- 
ally available ruby optical maser* which produces 
approximately 3 joules of 6943A light in a one- 
millisecond pulse. This light is passed through 
a red filter for the elimination of the xenon flash 
background and is then brought to a focus inside 
a crystalline quartz sample. The emergent beam 
is analyzed by a quartz prism spectrometer 
equipped with red insensitive Eastman Type 103 
spectrographic plates. A reproduction of the first 
plate in which there was an unambiguous indica- 
tion of second harmonic (3472 A) is shown in Fig. 
1. This plate was exposed to only one “shot” 
from the optical maser. We believe the following 
two facts, among others, rule out the possibility 
of artifact: 

(1) The light at 3472 A disappears when the 
quartz is removed or is replaced by glass. 

(2) The light at 3472 A exhibits the expected 
dependence on polarization and orientation sum- 
marized in Table I. 

Considerations of the photographic image densi- 
ty and the efficiency of the optical system lead 


us to believe that the order of 10" second har- 
monic photons were generated within the quartz 
sample per pulse. 

The production of a second harmonic should 
be observable in isotropic materials such as 
glass if a strong bias field were applied to the 
sample. This bias could be oscillatory, thus 
producing sidebands on the fundamental fre- 
quency and the harmonics. 

We would like to thank the staff of Trion Instru- 
ments, Inc., for their valuable and sustained co- 
operation in this work. 





*This work was supported in part by the U. S. Atomic 
Energy Commission. 
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*R. J. Collins et al., Phys. Rev. Letters 5, 303 
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3a, Javan, W. R. Bennet, and D. R. Herriott, Phys. 
Rev. Letters 6, 106 (1961). Even though the intensity 
of the gas device is very low compared with ruby masers, 
the gain in coherence volume and the potential improve- 
ment of focussing suggest that the gas maser may be 
comparable or even superior as a source for optical 
harmonics. 
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DIRECT DETERMINATION OF X-RAY REFLECTION PHASE RELATIONSHIPS 
THROUGH SIMULTANEOUS REFLECTION* 


M. Hart and A. R. Lang 


H. H. Wills Physics Laboratory, University of Bristol, Bristol, England 
(Received July 18, 1961) 


This Letter reports a successful direct experi- 
mental determination of x-ray reflection phase 
relationships through coherent dynamical inter- 
actions in simultaneous x-ray reflection. When 
a crystal is so oriented with respect to the inci- 
dent beam that both reflections h (=hkl) and h’ 
(=h’k'l’) are excited, the region in reciprocal 
space concerned is in the vicinity of the point 
of intersection of the three spheres Z,, Zp, and 
<p centered, respectively, on reciprocal lattice 
points O, h, andh’, the radii of the spheres be- 
ing k=”K where » is the crystal refractive index 
and K the free-space wave number for the radia- 
tion used. In this region the dispersion surface 
is made up of three sheets, and in the crystal 
there are nine wave vectors kp?) (g=O,h,h’; 
i=1,2,3). (For diffraction at low angles, polar- 
ization effects may be neglected; otherwise one 
would have to consider six sheets.) Clearly the 
situation is very complicated when all the dif- 
fracted waves are comparable in amplitude. How- 
ever, we have achieved our result in a straight- 
forward way by three steps. 

(1) High-resoluticn x-ray diffraction topo- 
graphs’? are used to provide adequate definition 
of angular and spatial relations. 

(2) A valuable simplification is gained by re- 
stricting attention to an investigation of the mod- 
ification of the dispersion surface of the h reflec- 
tion as one proceeds along the line of intersection 
of =, and 2» and approaches the intersection with 
Zp 

(3) For studying small modifications of the dis- 
persion surface we make use of the phenomenon 
of Pendellésung interference. It has been shown’’ 
that in x-ray projection topographs? of wedge- 
shaped crystals in the Laue case the Pendelld- 
sung fringe pattern is determined by the mini- 
mum diameter of the dispersion surface hyper - 
bolae, which is proportional to the structure am- 
plitude |F)!. The fringe spacing, for a given 
wedge angle, is inversely proportional to this 
diameter. In the present investigation, displace- 
ments of Pendellésung fringes on the topograph 
provide the measure of modifications of the dis- 
persion surface. 

The geometry of the experiments is explained 
in the stereographic projection, Fig. 1. The 
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FIG. 1. Diffraction geometry for simultaneous re- 
flection. 


crystal (not shown) is at the center of the sphere. 
Diffracted beams lie in the hemisphere above the 
paper, whereas the incident beam lies in the hem- 
isphere below the paper. Continuous and inter- 
rupted lines are used, respectively, on the hem- 
ispheres above and below the paper. The line AA’ 
represents the horizontal and the crystal is illu- 
minated by a narrow vertical ribbon of x rays. 
The angular range of the incident-ray bundle is 
represented by the interrupted rectangle, much 
exaggerated in the drawing. The x-ray source 

is effectively a point and the vertical range of 
incident rays is determined by the height of crys- 
tal illuminated. There is thus a one-to-one cor- 
respondence between position on the crystal and 
the vertical angle of the ray incident upon it. 

The plane normals h and h’ are shown, together 
with the cones generated by all rays making the 
appropriate Bragg angle with these planes for 
the radiation used (Ag Ka,). The plane normal 

h lies horizontal and the whole of the incident 
beam straddles the h cone. A diffracted beam 
the full height of the crystal is thus obtained, 
incident rays X and X’ giving rise to diffracted 
rays X and X’, respectively. Next, the crystal 
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is rotated about plane normal # to bring the h’ 
cone into the incident beam, but in this case the 
incident beam cuts only a small part of the cone. 
The incident ray common to both the h cone and 
h’ cone is at Z,, the small bundle of the h’ dif- 
fracted beam is at Z),, and in the extended h dif- 
fracted beam we find an Aufhellung® in the direc- 
tion Zy corresponding to the incident ray Z,. The 
crystal has the shape of a gently tapering wedge, 
with the thickness contours making a moderate 
angle with the horizontal. The projection topo- 
graph obtained by traversing the crystal back 

and forth appears as sketched in the left-hand 
part of Fig. 1. The thin end of the wedge is in 
the upper right and the Pendelldsung fringe order 
is in the sequence indicated. (The diffracted 
beam Z,, may possibly pass through the dif- 
fracted beam screening slits and appear as a 
streak as shown.) 

The electron diffraction dynamical theory of 
Kambe and Miyake,**” modified to the x-ray case, 
shows that in the vicinity of 2p), the minimum di- 
ameter of the / reflection dispersion surface hy- 
perbolae is proportional not to |F,| but to quan- 
tities IF ,| just inside 2p, and |F_| just outside 
Zp, to a good approximation. Whether |F,! is 
greater or less than |F _! depends upon whether 
cospp is positive or negative, Ppp’ being a 
function of the individual phase angles, a, of the 
reflections h, h’, andh’ -h, 


*. “5° 


Phy h'-h’ 


In a centrosymmetric cyrstal, pp can only be 
zero or 7; hence cos@p p,: is then restricted to 

the values +1. Now wave vectors of rays incident 
between X and Z, (and hence diffracted between X 
and Zy) lie outside =;-, whereas those incident 
between X’ and Z, lie inside Zp. Thus, on the 
topograph, one expects the effective diameter of 
dispersion surface hyperbolae ‘to be proportional 
to |F ,| just above the Aufhellung, and to IF _| 

just below the Aufhellung. Hence if |F | is great- 
er than |F _| the Pendellésung fringe spacing is 
decreased above the Aufhellung and increased be- 
low; and vice versa. In the former case the fringe 
discontinuity appears as sketched. The experi- 
ment is quite a sensitive one, since at the fringe 
of order 5, a change of only 10% in dispersion 
surface diameter will shift the fringes by half a 
period. 





FIG. 2. Diffraction topograph showing dislocation of 
Pendellésung fringes at an Aufhellung. 


A topograph corresponding to the sketch in Fig. 
1 is shown in Fig. 2. The crystal is germanium, 
the reflection h is 220 and the reflection h’ is 311. 
The Aufhellung appears across the field and the 
streaks due to the h’ reflection (both Ka, and Ka,) 
appear across the lower half. Although the crys- 
tal thi¢kness variation is not as smooth as could 
be desired, and numerous dislocations are scat- 
tered over the field, the form of the shift of the 
Pendellisung fringes clearly shows that |F | is 
greater than |F_]}. 

Since the structure is centrosymmetric, the 
observation shows that cosép,y,-=1, and hence 
Opp1=Ap - Apr t+Apr_p,=0. This result is in 
agreement with that calculated for the known 
structure. A direct experimental determination 
of phase relationships is thus effected. 

The authors would like to thank Professor F. C. 
Frank and Professor S. Takagi for advice. 
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European Office. 
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EXPERIMENTAL TEST OF THE SPIN-WAVE THEORY OF A FERROMAGNET 


A. C. Gossard, V. Jaccarino, and J. P. Remeika 


Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received July 24, 1961) 


A detailed theory of the thermodynamic proper - 
ties of an idealized model of a ferromagnet, in 
which there exist only temperature -independent 
ferromagnetic nearest neighbor interactions hav- 
ing the form J})/"8-§;, has been given.’ There 
has hitherto been no satisfactory test of the the- 
ory, primarily because no known ferromagnetic 
materials satisfied the assumptions of the model. 
Recently, Tsubokawa?” has discovered an insula- 
ting ferromagnet, CrBr,, which closely approx- 
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FIG, 1. Crystal structure of CrBrs; (from Tsubokawa). 
The Br” ions form a nearly perfect close-packed hexag- 
onal structure. The Cr** ions occupy two-thirds filled 
hexagonal layers with two Br~ layers separating each 
Cr* layer. Each Cr** is surrounded by a nearly per- 
fect octahedron of Br~ ions. At the right a simplified 
model used to calculate the spin-wave dispersion law 
is shown. 
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imates the model. Using nuclear magnetic reso- 
nance (NMR) techniques, we have made a precise 
measurement of the temperature dependence of 
the magnetization of this compound. 

CrBr, has the DO, (Bil,) hexagonal layer struc- 
ture® shown in Fig. 1 and is prepared by heating 
powdered chromium in an atmosphere of bromine 
at 800°C.? Tsubokawa has measured several of 
its magnetic properties. In particular he has 
found a ferromagnetic Curie point T; =37°K and 
a saturation magnetization of 3u B ber Cr** ion 
as would be expected for the 3d* *S configuration 
in an octahedral crystal field. For the purely 
ionic *S configuration a 375 000-oersted negative 
internal field at the Cr nucleus would be expected 
as a result of the exchange polarization of inner 
closed-shell s electrons by the aligned d elec- 
trons.‘ Since overlap of the Cr and Br wave 
functions and covalent bonding reduce the mag- 
nitude of the internal field at the Cr nucleus, we 
searched for the Cr** NMR in the frequency re- 
gion below 85 Mc/sec using a frequency-swept 
modified Robinson marginal oscillator.® 

An oscilloscope tracing of the spectrum ob- 
served at 58.038 Mc/sec (Hpsg = 241220 oe) and 
1.34°K in zero external field is shown in the up- 
per corner of Fig. 2. Three distinct absorptions 
are seen, whose separations and intensities cor- 
respond to the allowed transitions for J= 3 anda 
nuclear electric quadrupole interaction h*e?q,,Q 
= 1.184 Mc/sec. 

The occurrence of sharp satellites indicates a 
uniform orientation of the internal magnetic field 
with respect to the axes of the electric field gra- 
dient tensor. From this, and the external field 
dependence of the frequency for resonance, it 
may be concluded that the nuclei participating 
in the resonance are situated in the body of the 
domains and not in the domain walls. Enhance- 
ment of the resonant absorption occurs from 
domain rotation rather than wall motion.* The 
sign of the magnetic field at the nuclei was de- 
termined to be negative with respect to the direc- 
tion of magnetization from the observed decrease 
of the resonance frequency with increasing ex- 
ternal field. The temperature dependence of the 
Cr*®* NMR is shown in Fig. 2. The resonance 
was observed up to 15°K.” 

Because of the hfs interaction AI-5 between 
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the ionic spin $ and the nuclear spin i, NMR is 
observable at a temperature T<T,, at a frequen- 
cy v(T) =h*A(S)z, where (S)7 is the thermal 
average of §. Since M(T) =NgB(S)7, we expect 
v(T) to be proportional to M(T) provided that A 
depends neither implicitly nor explicitly on T. 
The temperature dependence of hfs interactions 
has been investigated in other magnetic solids.® 
The results of such studies lead us to conclude 
that the variation of A between 0° and 4°K, in 
CrBr,, is considerably less than the intrinsic 
precision of our experiment (1 part in 30000). 

The spin-wave prediction® for the temperature 
dependence of the spontaneous magnetization 
M(T), for T«T, is, 


M(T) =M(0)[1 -aT™ - oT” -cT” -dT* -eT™ ----]. 


(1) 


The half-odd-integral terms beyond T* result 
from the discreteness of the lattice, while the 

T* term is the lowest spin-wave interaction term. 
In the case of cubic symmetry with no external 
field or anisotropy, the coefficients a, b, c, d, 
and e are functions only of J, S and the lattice 
geometry. In the DO, lattice, however, several 
exchange interactions are involved. Nearest 
neighbor spins within a layer are strongly cou- 
pled by superexchange through a single inter- 


T 3/2 (° «3/2 ) 


vening ion. Exchange between spins in adjacent 
layers is weaker because two intervening ions 

in series separate the interacting spins. Calcu- 
lation of the expected magnetization coefficients 
is thus complicated by the fact that some weak 
exchange is expected even for spins not directly 
above one another. To simplify calculation of 
the coefficients, we have adopted a model of the 
lattice shown at the right in Fig. 1. Only two ex- 
change interactions, J; and Jj, are considered 
in the spin-wave calculation. For J7«J; we ob- 
tain a=0.0130 J,-7J,*and b =0.000140 JJ. 
In the presence of a magnetic or anisotropy field, 
H, the T* coefficient in the same formalism be- 
comes 


(a/2.612)2/ [n “#2 exp(-ngBH /kT) |. 
The other coefficients are modified similarly. 

In Fig. 2 the results of the NMR measurements 
are plotted against T** for T<0.117T,. At the 
lowest temperatures the effect of the anisotropy- 
modified T** term is seen to predominate. Using 
Dillon’s observed value’® of 6850 oersteds for 
the anisotropy field, a least-squares fit of the 
form a <.e MEBH/kT 

- ee. . Pet v2 
v(T) vo z61a~ a7 —7 
b 2 "ebH /kT 
-aq7 aT” ’ 
1.341 ‘ n 
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was made. v(0) was determined to be 58.096 Mc/ 
sec, corresponding to a field at the nucleus of 
241 460 oersteds. a was (2.544+ 0.067) x 1073 °*K~*2, 
b was (3.03 + 1.04) x10~°°K~¥?, and the mean square 
square deviation of the data points from the fit 

was 0.0016 Mc/sec. Inclusion of T” and T* 

terms does not significantly improve the fit. 

From a and 6 we calculate that J; =5.44°K and 
J,=0.88°K. The smallness of Jj relative to J;, 
inferred from inspection of the superexchange 
“paths,” is thus confirmed, and manifests itself 

in the magnetization curve as an unusually large 
T* term. 

We wish to thank R. J. Blume for the circuit of 
the Robinson-type spectrometer, J. L. Davis for 
experimental assistance, E. M. Kelley for aid in 
growing the crystals, and Dr. L. R. Walker for 
theoretical counsel. 
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PHASE TRANSITION IN MERCURY TELLURIDE™ 


J. Blair and A. C. Smith 
Department of Electrical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received July 17, 1961) 


An abrupt change in the resistivity of mercury 
telluride has been observed at a pressure of about 
15000 kg/cm? at room temperature. This change 
is ascribed to a transition between two polymor- 
phic phases of HgTe. The resistivity changes es- 
sentially reversibly by a factor of 10* to 10° at 
this transition. No information about the struc- 
ture of the high-pressure phase has been obtained. 

Figure ‘1 shows the variation of the resistivity 
of a sample of HgTe as a function of pressure. 

The resistivity was measured with a four-probe 
method using a precision potentiometer. At 16000 
kg/cm’, the resistivity increases by a factor 
greater than 10*, The transition takes place very 
slowly and the points above 16000 kg/cm? do not 
represent an equilibrium situation. The time that 
elapsed while the sample was in the high-pressure 
phase was several hours and the resistivity drifted 
upward consistently. As the pressure was de- 
creased, the reverse transition occurred at 12000 
kg/cm*, The resistivity decreased by a factor of 
10*, but did not quite return to its original value. 

When the sample was examined, several cracks 
were found. These cracks presumably account for 
the difference in the atmospheric pressure values 
of the resistivity before and after the application 
of pressure. The sample was initially a single 
crystal and if the transition nucleated at a num- 
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ber of points, the sample would be polycrystal- 
line after the transition. The strains arising from 
misorientation of the crystallites would produce 
the observed cracks. 

The resistivity of this sample before the nres- 
sure experiment was 1.5x10~* ohm cm. Meas- 
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urements by Nelson! on adjacent samples from 
the same ingot yield an acceptor concentration 

of 3.7x10'® cm~, and an electron mobility of 

2.3 x10* cm*/volt sec. The change in resistivity 
below the transition pressure is quite reproduc- 
ible and is essentially the same for samples from 
different ingots. A detailed explanation of this 
behavior is planned for a later paper. 

The phase transition was observed in another 
sample of HgTe and an effort was made to de- 
termine the transition pressure more accurately. 
The limits of the transition pressure were again 
found to be 16000 kg/cm? with increasing pres- 
sure and 12000 kg/cm? with decreasing pressure. 


The authors wish to express their thanks to 
Professor William Paul and the High Pressure 
Group at Harvard University for the use of the 
pressure equipment. We also thank Dr. R. E. 
Nelson for supplying the HgTe samples. 





*This work was partially supported by the Air Force 
Cambridge Research Laboratories (Air Research and 
Development Command), U. S. Air Force. 
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NEGATIVE PION-PROTON ELASTIC SCATTERING AT 1.51, 2.01, AND 2.53 Bev/c 
OUTSIDE THE DIFFRACTION PEAK* 


Kwan W. Lai, Lawrence W. Jones, and Martin L. Perl 


The University of Michigan, Ann Arbor, Michigan 
(Received July 27, 1961) 


The differential elastic scattering cross sec- 
tions for negative pions on protons have been 
measured at incident pion momenta of 1.51, 2.01, 
and 2.53 Bev/c with emphasis on the angular re- 
gion outside the diffraction peak. The purpose 
of the experiment was to examine the behavior 
of the large angle differential elastic cross sec- 
tion as a function of energy from the energy of 
the highest known resonance in the pion-nucleon 
system into the region where the total cross sec- 
tions appear to be approaching an asymptotic val- 
ue.? 

The experiment was performed at the Bevatron 
using a luminescent chamber system to photo- 
graph the tracks of the scattered pion and the 
recoil proton from a liquid hydrogen target. The 
chamber configuration used is shown in Fig. 1, 
where the plane of the figure is normal to the op- 
tic axis of the image intensifier tube system. 
Lenses focused this image of the sodium iodide 
scintillators onto the first image tube cathode, 
and the gain of the system was sufficient to re- 
cord single photoelectrons from this cathode. 

A coincidence between the coplanarity counters, 
which subtended Aé =+7.15° from the vertical 
plane, triggered the system on probable elastic 
events. This post-event triggering enabled the 
tracks to be subsequently recorded on film with 

a time resolution of a few microseconds. A de- 
scription of the apparatus has been reported else- 
where in detail.? The film was double scanned, 


and angles for elastic scattering events were re~- 
corded with suitable corrections for distortions 
in the image tube system. The angles and ioni- 
zation of the two tracks were used to identify an 
event as elastic, and subsequent correlations of 
the track angles with kinematics were consistent 
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FIG, 1. Geometry of liquid hydrogen target, scintil- 
lation counters, and sodium iodide luminescent cham- 
bers. The coplanarity scintillation counters C; and C) 
subtended +7.15° from the vertical plane (the plane of 
the drawing) for scatterings on the beam axis. Besides 
anticoincidence counters A; and A», two other anticoin- 
cidence counters flanked the target parallel to the plane 
of the figure to bias against inelastic events. S3 is the 
final beam defining counter. The sodium iodide crys- 
tals were imaged by lenses onto the image tube system 
which was oriented perpendicularly to the plane of the 
figure. 
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FIG, 2. Elastic differential cross section of ™~+p 
~~ +p in the center-of-momentum system. 6* is the 
angle of the scattered pion in that system, The values 
of d0/dQ* at the smallest angle measured at each ener- 
gy are too large to be placed at the correct places on 
the graphs, but their values are indicated. 


with the known angular resolution of the measure- 
ments. 
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2405 elastic scatterings were analyzed at 1.51 
Bev/c, 1300 events at 2.01 Bev/c, and 1080 
events at 2.53 Bev/c. The results are presented 
in Fig. 2. The differential cross sections have 
been corrected for multiple Coulomb scattering 
and nuclear interactions in sodium iodide, for 
geometrical and scanning biases, and for pseudo- 
elastic events from the target vessel (from data 
taken with the target empty). The absolute value 
of the differential cross section may be uncertain 
by as much as +15% from systematic errors as- 
sociated with this new technique. The data at 1.51 
Bev/c are in good agreement with the bubble 
chamber result of Chretien et al.° 

From the data it may be noted that the back- 
ward hump, which has a maximum height of 2.1 
mb/sr at 900 Mev and 1.1 mb/sr at 1020 Mev,‘ 
is down to 0.4 mb/sr at 1.51 Bev/c (1.37 Bev) 
and is not present at 2.01 or 2.53 Bev/c. The 
angular distributions behind the diffraction peak 
at 2.01 and 2.53 Bev/c are roughly constant with 
angle, decreasing from 0.18 mb/sr at 2.01 Bev/c 
to 0.11 mb/sr at 2.53 Bev/c. Although the data 
at these two higher momenta can be taken to sug- 
gest some oscillatory structure behind the dif- 
fraction peak, they are not inconsistent with an 
isotropic distribution behind the diffraction peak. 

We wish to acknowledge our appreciation to 
Dr. E. J. Lofgren, his associates, and the staff 
of the Bevatron of the Lawrence Radiation Labo- 
ratory for their hospitality and assistance through 
out the preparation, execution, and analysis of 
this experiment. 
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TOTAL CROSS SECTIONS FOR PIONS ON PROTONS IN THE MOMENTUM RANGE 4.5 TO 10 Gev/c 


G. von Dardel, R. Mermod, P. A. Piroué,* M. Vivargent, G. Weber, and K. Winter 


CERN, Geneva, Switzerland 
(Received July 17, 1961) 


In a previous paper’ we reported on measure- 
ments of total cross sections for p, 6, K*, and 
a on protons at high energy. According to a 
theorem given by Pomeranchuk,?’ the total cross 
sections for a particle and its antiparticle must 
become equal if they approach constant values 
at high energy. For pions, the conditions of this 
theorem are thought to be fulfilled at lower ener- 
gies than for the other particles. In our previous 
paper, a large uncertainty on the muon contami- 
nation did not allow us to draw any significant 
conclusion concerning the validity of the Pomer- 
anchuk theorem. 

The present paper reports on a more accurate 
measurement of the (7+, p) total cross sections 
in the momentum range 4.5 to 10 Gev/c; the ex- 
perimental method used eliminated the effect of 
muon contamination on the cross sections. 


The layout of the experiment is shown in Fig. 1. 


A beam of secondary particles enters the experi- 
mental hall and comes to a focus at the center of 
a velocity-selective gas Cerenkov counter C, af - 
ter having traversed a system of two quadrupole 
lenses and a bending magnet located in the target 
area. The beam spot at the Cerenkov counter is 
imaged by another doublet of quadrupole lenses, 
Q;,, Q,, onto transmission counters, S,, S,, S,, 
which are located behind a 312-cm long liquid 
hydrogen target. Bending magnet B, provides a 
second momentum analysis eliminating secon- 
daries from reactions produced in the collimator 
walls, in the Cerenkov counter, and in the air, 
and also eliminating part of the decay muons. 
The beam incident on the hydrogen target is de- 
fined by scintillation counters, S,, S,, and S,, 2, 
5, and 2 cm in diameter, respectively, operated 
in coincidence with the Cerenkov counter, which 
is set to count pions, but does not discriminate 
against muons and electrons. After traversing 
the transmission counters S,, S,, and S,, having 
diameters of 10, 15, and 20 cm, respectively, 
the beam is intercepted by a 1-mm long iron 
absorber in which strongly interacting particles 
and electrons are lost, whereas more than 99% 
of the muons are transmitted into counter S, (30 
cm in diameter). The pions transmitted into 
counters S,, S,, and S, are signalled by the co- 
incidences (S,S,S,CS,5,), (S,S,S,CS,5,), and 
(S,S,8,CS,5,), respectively. These rates, nor- 
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FIG, 1. Experimental layout. The first doublet of 
quadrupoles and the first bending magnet are in the tar- 
get area (not shown). Lead collimators shown are 2.5 
x2.5 em’, S,-S, scintillators, Q34 quadrupole lenses, 
B, second bending magnet. 


malized to the incident rate (S,S,S,C), were meas- 
ured for both target “full” and target “empty.” 
Two steel plates of appropriate thickness simu- 
lated the empty target. The coincidence logics 
was arranged such that the time resolution in 

all channels was determined by threefold coin- 
cidence units of 5-nsec resolution. The acciden- 
tal rates were negligibly small. 

The electron contamination in the negative beam 
was estimated to 0.6% at 6 Gev/c from a meas- 
urement using a threshold Cerenkov counter in 
a similar beam, and a calculation of the effect 
of radiation losses in the material in the actual 
beam used in the experiment. All sources of 
positron contamination are expected to yield a 
smaller fraction in the positive beam. In both 
beams, the contamination is decreasing with in- 
creasing energy. The correction to the cross 
section and to the difference o(1~, p) - o(n*, p) 
will be less than 0.2 mb. 

Earlier measurements of the transmission in 
hydrogen of negative pions had shown that there 
is no structure within 3% of the cross section 
for either isospin state in the momentum range 
5-10 Gev/c. Therefore, it was sufficient in the 
present work to measure cross sections at a few 
momenta only (4.5, 5.75, 7.0, and 10.0 Gev/c). 

Total cross sections were obtained by extrap- 
olating to zero solid angle the cross sections 
measured with counters S,, S,, and S,, subtend- 
ing solid angles averaged over the target length 
of 0.469, 0.826, and 1.315 msr, respectively. 
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Table I. Total cross sections for negative and posi- 
tive pions on protons and their difference. 











Momentum  o(1~,p) o(n* ,p) o(1~ ,p) -a(n* ,p) 
(Gev/c) (mb) (mb) (mb) 
4.5 30.2+0.4 27.6+0.3 2.620.5 
5.75 29.120.4 26.9+0.4 2.2+0.6 
7.0 28.4+0.4 26.1+0.3 2.3+0.5 
10.0 27.0+0.4 25.2+0.4 1.8 +0.6 








For the lowest momentum, a correction for mul- 
tiple scattering of 1.3% was applied to o(S,). The 
cross sections plotted versus solid angle form a 
straight line at all momenta, and consequently, 
linear extrapolation to zero solid angle was used 
to derive the total cross section. The results are 
given in Table I. 

The errors given are statistical only, including 
the error in extrapolation. The absolute values 
of the cross sections may have an additional er- 
ror of about 0.5 mb due to uncertainties in the 
effective target length and the calibration of the 
steel plates simulating the empty target. This 
should affect neither the difference of cross sec- 
tions for negative and positive pions, nor the mo- 
mentum dependence. 

The results are shown in Fig. 2 together with 
data of other authors.*” Our value for the (1*, p) 
cross section at 4.5 Gev/c is lower than values 
given by Longo et al.* and Chzan et al.” 

Unlike the situation for (p, p) total cross sec- 
tions® which have been shown to be constant be- 
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FIG, 2. Total cross sections for (™-,p) and (7* ,p) 
versus momentum. Observe break in ordinate scale. 


128 


tween 5 and 28 Gev/c, a monotonic decrease is 

observed for both the (7*, p) and (1~, p) cross sec- 
tion with increasing energy. Hence, the condition 
for Pomeranchuk’s theorem? that total cross sec- 
tions tend asymptotically to constant values is not 
satisfied for the 7-nucleon system up to 10 Gev/c, 

The difference observed between the cross sec- 
tions for negative and positive pions is outside the 
errors, and slowly decreasing with increasing en- 
ergy. This difference may be due to elastic charge 
exchange of negative pions, but should vanish if the 
equality of the cross sections in different isotopic 
spin states, predicted by Pomeranchuk’s theorem, 
were fulfilled. That it does not vanish confirms 
the conclusion of the preceding paragraph. 

The monotonic decrease of the cross sections at 
high energies could, if it continued, be considered 
in connection with recent theories by Berestetsky 
and Pomeranchuk’® and Gribov!® which suggest a 
cross section decreasing to zero faster than 
(InE)~, Extrapolating our values according to 
the functional dependence o=o,[In(E/E,)]~ leads 
to cross sections of 20 mb at 100 Gev and 16 mb 
at 1000 Gev. 

We are indebted to the Machine Group of the 
CERN Proton Synchrotron for the efficient opera- 
tion of the machine during the whole run. We also 
thank M. Rousseau and L. Velati for their assist- 
ance in setting up and running the experiment. 
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MAGNETIC MOMENT OF POSITIVE AND NEGATIVE MUONS* 


D. P. Hutchinson, J. Menes,f and G. Shapiro 
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A. M. Patlach 
International Business Machines Watson Scientific Computing Laboratory, New York, New York 


and 


S. Penmant 


University of Chicago, Chicago, Illinois 
(Received July 24, 1961) 


This Letter reports on the refinement of the 
technique for measuring muon magnetic moments 
and on its application to positive and (bound) neg- 
ative muons in various materials. These refine- 
ments enable us to obtain a statistical accuracy 
of several parts per million in reasonably short 
runs. In addition to the absolute measurement of 
Su» this sensitivity permits the observation of 
environmental effects and thus has applications 
to nuclear and solid-state physics. 

The magnetic moment of the u» meson, like that 
of the electron, is of theoretical importance be- 
cause of the exact predictions that can be made 
of its magnitude. Quantum electrodynamics en- 
ables one to calculate the moment? in units of 
the muon magneton, ef /m uc: TCP invariance 
requires the positive and negative muon, as par- 
ticle and antiparticle, to have equal magnetic 
moments.” 

Since the discovery of parity nonconservation, 
there has been a succession of increasingly ac- 
curate measurements*® of the muon magnetic 
moment, at this laboratory and elsewhere. These 
experiments, as well as the present one, gener- 


_ally express the muon moment in units of the pro- 


ton magnetic moment, and depend on separate de- 
terminations of the muon mass*® for comparison 
with theory. A direct measurement of the “anom- 
alous part” of the magnetic moment, which does 
not depend strongly on the mass value for its in- 
terpretation, has also recently been made.” The 
accuracy of the present experiment so far ex- 
ceeds that of any muon mass determination that 

it must await future experiments on the latter 


before its significance can be realized. Alterna- 
tively one can accept present theoretical’ or ex- 
perimental’ estimates of the g factor, and com- 
bine them with the present results to obtain a 
“best” value for the muon mass. 

The basic method of this experiment is similar 
to that of reference 5, to which the reader is re- 
ferred for details. The increased sensitivity of 
the present experiment is due to numerous im- 
provements in the apparatus, including a higher 
muon precession frequency, an extremely homo- 
geneous magnetic field, and the transistorizing 
of all circuitry. Recent improvements in the 
beam facilities at Nevis, including the long-duty- 
cycle beam and a new low-energy muon beam, 
also contributed to the increased sensitivity. Fig- 
ure 1 shows the experimental arrangement. 

Figure 2 shows the results of typical runs with 
positive and negative muons stopping in HCl solu- 
tion and graphite, respectively. Early-late phase 
shift is plotted against magnetic field, as meas- 
ured by a proton magnetic resonanc® at the moni- 
tor position. Experimental conditions were such 
that for a target material that does not depolarize 
positive muons (e.g., bromoform), running time 
of three hours was sufficient to reduce the sta- 
tistical error to 12 ppm. 

Results of the runs with positive muons stopping 
in various target materials are summarized in 
Table I, in chronological order. Statistical er- 
rors calculated from a priori considerations are 
indicated. 

The accuracy of the experiment is such that di- 


amagnetic shielding of the fini, by djomic elec- 
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Table I. Results for positive muons. 


Errors indicated are statistical only. 





Stopping material 


Proton frequency at 
resonance (kc/sec) 





Methylene iodide—CHyl, 
Distilled water —H,O 

HCl 

H,O-—repeat 

Methylene iodide — repeat 
Bromoform 


Average 
Monitor-to-center field difference 


Muon precession frequency at resonance 
/ 
I 'p 


»* in graphite (99.995 % pure electromolded) 


55 824.45 +0.49 
55 823.41 +0.49 
55 822.07 +0.35 
55 821.39 +0.94 
55 823.33 +0.59 
55 823.30 +0.69 
55 823.01 +0.43 

4.85 +0.10 
55 818.16 


177688.4 +0.1 
3.18334 +0.00005 


55 801.8 +1.8 





trons must be considered. If the positive muon 
substitutes immediately for a proton and experi- 
ences the same shielding, the shifts are expected 
to be less than 3 ppm (relative to a proton in wa- 
ter) for all liquid targets used.® The proton reso- 
nance for measuring the field is also taken in 
aqueous solution, so essentially no net correction 
has to be made for the ratio of the moments. 

It is noted that the substituted hydrocarbons 
yield consistently higher values than the aqueous 
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targets, beyond the predictions of chemical shifts. 
Moreover, there is a downward trend with time 
for runs repeated with the same type of target. 
Rigorous checks of the system gave no evidence 
for a systematic drift in this direction. There- 
fore, it is assumed that both these effects are 
statistical fluctuations and the error is increased 
to account for the observed spread. 

The average proton resonance frequency thus 
obtained is 55823.01+ 0.43 kc/sec. From this 
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FIG, 2. Early-late phase shift vs proton resonance 
frequency at monitor position. The straight line is the 
least-squares best fit, crossing the x axis at the field 
setting where the muon precession frequency equals 
the reference value. The horizontal flag at this point 
is a calculated a priori statistical error. (a) Positive 
muons in HCl solution; (b) negative muons in graphite 
(reactor grade). 


must be subtracted 4.854 0.10 kc/sec as the dif- 
ference between the field measured at the moni- 
tor position and at the center of the target. The 
muon precession frequency at “resonance” is 
177 688.4+0.1 kc/sec. An improved transistor 
regulator’® permits maximum excursions of less 
than 5 ppm in the field setting (~ 13.4 kgauss) dur- 
ing each one-hour run. Field spatial homogene- 
ity is such that the maximum deviation from the 
average over 95% of the target volume (2 in. x3 
in. X1 in.) is not more than 10 ppm. The result 
for positive muons is f, If = 3.18334 + 0.00005, 
where the error has been increased to allow for 
possible systematic errors. 

Similar measurements have been made on the 
magnetic moment of the electron. The most ac- 
curate experimental result to date is** Se/tp 
= 658.22759 + 0.00004, as the ratio of the spin 
precession frequency of free electrons to that 
of protons in water. Since in the present experi- 
ment both muon and proton are bound, a correc- 
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tion’? for diamagnetic shielding of 25.6 ppm must 
be made to one of these results before they can 
be compared. Making it to the electron experi- 
ment yields fe/fp(both free) = 658.2107, with an 
error small compared to that of the muon result. 
These frequencies are proportional to the g fac- 
tors of the particles, and inversely to their mas- 
ses. For the ratio of the g factors, one can use 
either the theoretical value, g, /gg = 1.000006; . 
or the experimentally determined g factors, 

8p =1.001145+ 0.000022,” and ge = 1.0011609 

+ 0.0000024.'* Using the theoretical g factor, 

one obtains for the ratio of the masses m,,/me 

= 206.768+ 0.003. Using the experimental num- 
bers, one has m,,/me =206.763+0.005. The 
larger error in the latter case reflects the ef- 
fect of folding in the error in the result of Char- 
pak et al. with that of this experiment. The av- 
erage of the mesonic x-ray determinations of 

the mass® yields m,,/m, =206.76+ 0.02, in excel- 
lent agreement with the above, but not sufficiently 
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Table II. Results for negative muons. Errors indicated are statistical only. 
Theoretical corrections (*10*) 
Binding 
Proton frequency effect on Nuclear Diamag- 
Stopping at resonance (< ~£ } x 108 Direct, radiative _polari- netic Knight 

Z material (ke/sec) a binding correction zation’ shielding shift’ Total 
u ” (corrected for 

diamagnetism) 55 821.5 +0.5 
6 Graphite, 

reactor grade 55 863.4 +1.7 - 7.540.3 - 6.29 -0.08 +0.04 -1.99 ? = 8.32 
6 Graphite, 99.995 % 

pure, 300°K 55 860.5 +3.1 - 6.9+0.6 - 6.29 -0.08 +0.04 -1.99 ? = 8.32 
6 Graphite, 99.995 % 

pure, 77K 55 865.7 +2.7 - 7.9+0.5 - 6.29 -0.08 +0.04 -1.99 ? = 8.32 
8 Oxygen in 

water 55 873.7+5.7 - 9.3+1.0 -11.04 -0.13 +0.12 -3.25 0 -14.30 
12 Magnesium 55 968 .5+3.7 -26.320.7 -23.79 -0.29 +0.53 -6.29 ? =-29.84 
14 Silicon 300 ohm cm, 

n type 56 023.5 +6.0 -36.121.1 -31.72 -0.40 +0.90 -7.95 ? -39.17 
16 Sulfur 56 090.0+9.1 -48.1+1.6 -40.35 -0.51 +1.4 -9.70 0 -49.16 








a 
See reference 17. 
See reference 19. 
See reference 15. 


accurate to afford a real comparison. 

A run was made with positive muons stopping 
in graphite. As noted in Table I, there is a para- 
magnetic shift of about 400 ppm from the average 
of the liquid targets. A similar (but smaller) 
shift was earlier™ noted in copper and aluminum, 
and was then interpreted as a Knight shift.** No 
other explanation of this effect is known to the 
authors; further investigations are being made. 
The graphite result was not used in the average 
for the positive muons. 

The results for negative muons stopping in 
various target materials are tabulated in Table 
II. These muons are observed in the ground 
state of mesonic atoms formed with the nuclei 
listed, and corrections must be made in order 
to compare with positive muons. The largest 
correction is due to binding.*® The statistical 
error on the data requires that the finite nuclear 
radius be taken into account in calculating this 
correction for Z 212. This experiment thus of- 
fers another way of probing the nucleus. The 
numbers quoted in Table II are from Ford, 
Hughes, and Wills.*” 

The next largest correction to be applied to 
negative muons is due to diamagnetic shielding 
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at the nucleus by the bound electrons.” Because 
of uncertainty in the final electronic states, Dick- 
inson’s’* results for atomic number Z-1 have been 
quoted without modification, although his results 
are for neutral isolated atoms. The magnitude 

of the effect is sufficient for this experimental 
method to provide a determination of the absolute 
shielding. 

Several smaller corrections taken from Ford, 
Hughes, and Wills’’ are also listed in Table II. 
No entries have been made under “Knight shift” 
for carbon, magnesium, and silicon. We hope 
to have soon a measurement of Knight shift in 
our carbon samples. The Knight shift in mag- 
nesium, expected to be a few parts in 10°, would 
predict a result in great disagreement with the 
data. A search was made for a “resonance” in 
the predicted region, with negative results. 
There exists a result for the Knight shift in 
silicon.”° However, the shift may be sensitive 
to the doping of the sample, and therefore this 
result has not been used. 

The negative-muon results are less accurate 
than these for positive muons because of depolar- 
ization and, for higher Z, shorter lifetime. 

Table II shows reasonable agreement between 
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predicted and measured muon moments in all the 
atoms investigated. It is concluded that, allowing 
also for theoretical uncertainties, the magnetic 
moment of positive and negative free muons are 
equal to within five parts in 10*. It may be that 
the various shifts mentioned will become a sub- 
ject for investigation in their own right, the mu- 
on serving as a tool for probing them. 
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Precise measurements of the magnetic mo- 
ments of negative 1 mesons bound in u -mesonic 
atoms are reported in the preceding Letter. 
Theoretical calculations of the magnetic mo- 
ments of u.-mesonic atoms are outlined and sum- 
marized in this Letter. 

It has been pointed out that the magnetic mo- 
ment (or g value) of a negative muon bound in an 
atomic orbit about a nucleus will be slightly less 
than that of a free muon and that the magnitude 
of this alteration is sensitive to the size of the 
nucleus.’ The corrections to the gyromagnetic 
ratio of a Dirac particle (g,=2) bound in the 
field of a zero-spin nucleus and surrounded by 
a zero-spin electronic cloud may be listed as 
follows: 

1. radiative correction, g,; 

2. binding correction to radiative correction, 


3. direct binding correction, g,; 

4. nuclear polarization correction, g,; 

5. electronic polarization correction, g,; 

6. electronic diamagnetic shielding correction, 
&e3 

7. center-of-mass correction, g,. 
We consider that the muon is bound in a 1s orbit 
about the nucleus. 

The radiative correction’ and the binding cor - 
rection® thereto have already been calculated: 


£,/Z,=a/2n +0.75(a/m)? = 0.001165, (1) 
&2/B = (26/150) a V) /mc’, (2) 


where (V) is the expectation value of the poten- 
tial energy and m is the muon mass. The direct 
binding correction, g,, is computed from a rela- 
tivistic treatment of the magnetic moment of a 
Dirac particle in a central field.‘ The result is 


83/8 = -(4/3) [Fdr, (3) 


where F is the small component of the radial 
wave function of the muon in the Coulomb field 
of a nucleus with finite size. Both the correc- 
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tions g, and g, have been obtained from the rela- 
tivistic bound-state muon wave functions calcu- 
lated on an IBM-704 computer. The nuclear 
charge distribution is taken to be 


p(r) = (Ze/Anr2N,) - aie for x <1, 
po 


=(Ze/4nr,°N,) (te for x >1, (4) 


where x =r/r, and 
N,=3+2n 7 +ne™, 


The parameters 7, andm, which characterize the 
nuclear radius and surface thickness, were cho- 
sen to be those which have been used to fit high- 
energy electron scattering results.* The effect 
of the nuclear size on g, is significantly larger 
than the experimental error for Z 212. 

A significant nuclear polarization correction 
term, g,, arises in second order perturbation 
theory from the magnetic interaction between 
the muon and the nucleus together with the mag- 
netic interaction of the nucleus in the external 
static magnetic field. In an approximate calcu- 
lation the term g, can be related to M1 nuclear 
gamma transition rates. An electronic polariza- 
tion correction term, g,, arises in second order 
perturbation theory from the magnetic interaction 
between the muon and the electrons together with 
the magnetic interaction of the electrons in the 
external static magnetic field.° The correction 
term g, is negligible unless the ground electronic 
state has fine-structure levels. 

The usual diamagnetic shielding’ by the atomic 
electrons will be very nearly the same for the 
muon as it is for the nucleus. This effect can be 
expressed by a correction term, g,: 


86/8 = “3a°Q) (ao/r,), 


in which the sum is taken over all the electrons. 
The center-of-mass correction g, is negligibly 
small. 
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Table I, Calculated corrections to the g values of selected u-mesonic atoms. 
a b 
Element 81/8 82/80 83/80 &/80 86/80 8/8 (g-gy+)/80 
gC 0.001165 -0.000008 -0.000629(0) +0.000004 -0.00019 1+0.00034 -0,00082 
30 0.001165 -0.000013 -0.001104(1) +0.000012 -0.00032 1-0,00026 -0.00143 
12Mg 0.001165 -0.000029 -0.002379(6) +0.000053 -0.00062 1-0,00181 -0.00298 
148i 0.001165 -0.000040 -0.003172(10) +0.000090 -0.00079 1-0.00275 -0.00391 
16S 0.001165 -0.000051 -0.004035(15) +0.00014 -0.00096 1-0.00374 -0.00491 








a s alas cerca , , 
Numbers in parentheses are uncertainties arising from 2% uncertainty in nuclear radius. 


The quantity g =i = 08; ° 


Table I lists the significant corrections for the 
elements measured experimentally. The terms 
g, and g, depend on the electronic state of the y- 
mesonic atom. For the elements listed, binding- 
energy considerations suggest that the » -mesonic 
atoms will be positive ions with Z-2 electrons and 
the electronic states will be 'S, except for oxygen 
and sulfur which will be *P,. The diamagnetic 
shielding term g, is obtained from Dickinson’s 
value® for an atom with Z-1 electrons together 
with a small estimated correction to give g, for 
the positive ion. The term g, will be negligible 
except for oxygen and sulfur where it amounts 
to about 5% and 2%, respectively; it is not in- 
cluded in the table since the large theoretical 
values for O and S are clearly inconsistent with 
the experimental values, thus indicating that the 
atomic states are not simple *P, states. 

In view of uncertainties about the electronic 
state of the u-mesonic atom (the state might 
even have electronic angular momentum J non- 
zero and then the hfs interaction between muon 
and electrons would be important), about chem- 
ical effects,° and about solid-state effects such 
as the Knight shift,’° the agreement between the- 
ory and experiment is close. The chemical shift 
and the Knight shift are paramagnetic effects, 
having the same sign as that required to account 


for the small remaining discrepancies. There 
is no indication for any anomalous behavior of 
the muon such as its polarization in the strong 
electric field of the nucleus. 

A more detailed report on these theoretical 
calculations will be submitted to The Physical 
Review. 
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EXAMPLE OF THE DECAY A°+p+p +0* 
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The decay A°+p+u~+¥, while expected, has 
not yet been reported to the knowledge of the au- 
thors. We report here a clear case of such a de- 
cay found in an exposure taken in the 20-in. BNL 
hydrogen bubble chamber. 

The exposure was taken in a carefully con- 
structed 7~ beam of energy 900 Mev and resolu- 
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FIG, 1. (a) Photograph of the event. (b) Tracing of 
the pertinent tracks (reduced). 


136 


tion 0.15%. The beam energy has been further 
confirmed by the analysis of over 1000 normal 
A°K® events. A picture of the event is shown in 
Fig. 1. The incident beam track, (1), interacts 
at a position 13 cm inside the chamber entrance 
which corresponds to an energy of 900 Mev. The 
A°, (2), goes 5.6 cm before decaying, a distance 
less than 2 mean lives. Track (3) is a proton. 
The w~, track (4), goes 9.6 cm before coming 
to rest and decaying into an electron, (5), which 
leaves the chamber through the top window. 

Four views of the event are available and were 
used in the measurement. The measurements 
and assigned errors are given in Table I. 

From the production angle of 16°+ 0.8° and 
known beam energy, the A° momentum is 495735 
Mev/c. Under the assumption of the decay A° 
>p+u~+v, the calculated neutrino energy from 
momentum balance is 79+ 15 Mev/c, giving a 
total final A° energy of 1223 Mev+ 20 Mev from 
summing the energy of its decay products. This 
energy and the momentum taken above, yield a 
mass for the A° of 1118+21 Mev, in good agree- 
ment with the known A° mass. 

The competing decay sequence A°+7~+p, 7” 
+u~+D is found to be very unlikely. The major 
points against this interpretation are as follows: 

(a) The proton momentum would have to be 475 
+20 Mev/c for this interpretation to hold. It is 
measured as 412+ 20 Mev/c. 

(b) If one assumes a pion to have been pro- 
duced, then its direction would have been such 
as to make too large an angle with the yu, by 
about twice the measurement error. 

(c) There is less than 1 mm of possible flight 
path for the 7~ since no sign of it is detected 
visually. The mean decay path for a 7~ with 
the required energy would be 5 meters. 

The above combine to give this interpretation 
a probability of occurrence of the order of less 
than 1075, 

Other possibilities for this event are: 


(1) It originates at the wall of the chamber, 
rather than at the observed pion ending. It would 
have had to travel 5 mean lives before decay. 

(2) The A° scatters in the chamber before de- 
cay, with no visible recoil. 

Each of these possibilities seems even more 
remote than the one previously explored since 
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Table I. Summary of measured data. 








Tracks Length measured (cm) Azimuth (deg) Dip (deg) Momentum (Mev/c) 
D5, (1) 2.5 62.5 +0.5 89.8 +1.0 1030 +2 
A° (2) 5.6 50.6 +0.2 100.8+0.4 495t}5* 
Pp (3) 23.0 56.3 +0.3 110 +0.4 412 +20 
(4) 9.4 50.8 +1.5 40.5+1.5 56 +1 (range) 





* Fitted quantity. 


one must still require the decay of the pion with- 
in 1 mm. 

From these different results, we conclude that 
this event has to be interpreted as a muonic de- 
cay of the A° particle. Theoretically, the rate 
for this decay mode is predicted to be 2.4x107°.! 

The above-described event is one of a sample 
of 3000 A°’s presently being analyzed. Our sys- 
tematic study of these events has as yet not been 
completed and we cannot set an upper limit on 
the rate of the decay. 

We would like to express our gratitude to Dr. J. 
Steinberger for valuable discussions. We would 


also like to thank Dr. R. Shutt and the members 
of the BNL Bubble Chamber Group for their in- 
valuable aid in obtaining the pictures and Mrs. 
Frances Joyce for having found the event. 
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CHARGE ASYMMETRIES IN THE ANGULAR DISTRIBUTION OF 7 AND K MESONS 
FROM ANTIPROTON ANNIHILATIONS IN FLIGHT 


B. C. Maglic, G. R. Kalbfleisch, and M. L. Stevenson 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received July 27, 1961) 


Presented here are the results of our search 
for a difference in the angular distributions of 
negative and positive 7 and K mesons produced 
in pp annihilations, at the center-of-mass mo- 
mentum of 657+16 Mev/c (1.6 Bev/c in the lab). 
There have been no direct theoretical predictions 
for such a difference. The statistical model pre- 
scribes an isotropic distribution for mesons of 
all charges.’ However, the model of Koba and 
Takeda? assumes that 2 or 3 annihilation pions 
are emitted from the antiproton and proton clouds. 
For annihilations in flight, the instantaneous mo- 
mentum of the pions from the f cloud will have a 
component parallel to the momentum, in addi- 
tion to their internal momenta. Since the model 
assumes that the pion cloud does not interact with 
the annihilation of the cores, it is reasonable to 
expect this component of the pion momentum to 
produce an excess of 7~ over 7* in the antiproton 
direction (forward). An equal excess of 7* should 


also be expected in the direction of the proton 
(backward). Evidence for such an effect has been 
scarce prior to this experiment. Angular distri- 
butions of pions of both charges from pp annihila- 
tions in propane at a c.m. momentum of 470 Mev/c 
(1.05 Bev/c lab) have not shown statistically sig- 
nificant deviation from isotropy.* A study of pp 
annihilations into two pions* at our momentum is 
also consistent with an isotropic distribution (8 7~ 
emitted forward, and 12 7~ emitted backward). 
However, annihilations into K*K™ pairs* suggest 
forward peaking for K~ (8 of the 11 K~ were emit- 
ted in the forward direction).° 

Our analysis was done on two samples of events: 
(i) 1620 events of the type b+p-+nz7, with n=4, 5, 
..., 8; and (ii) 287 events of the type j+p+K+K 
+n7, withn=1, 2, 3, and 4—in which at least one 
of the K mesons was observed to decay in the 
chamber. These events were observed inside a 
central volume of the 72-in. hydrogen bubble cham- 
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ber. The average antiproton momentum at the 
point of interaction was 1.61+0.04 Bev/c (925- 
Mev kinetic energy) for Sample (i), and for 86% 
of Sample (ii). The remaining 14% of Sample (ii) 
was from our run at 1.99+0.05 Bev/c. The ex- 
perimental conditions have been described previ- 
ously.® 

(i) Reactions )+p+nz7. Only annihilations into 
4 prongs have been analyzed, in which the identi- 
fiable reactions are 





p+p + 277+ 2n*+7°, (1) 
and 
p+p+2n7 +2n", (2) 


The abundances of these reactions are about 50% 
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FIG. 1. Histograms of c.m. angular distribution of 
pions from Reactions (1) to (3). (a) and (b): 7 and 1* 
from Reaction (1); (c): combined distributions (a) plus 
reflected (b); (d): 7~ and reflected 7 from Reaction 
(2); (e): negative and reflected positive prongs which 
did not fit Reactions (1) and (2); and (f): combined 
data (c), (d), and (e). The forward-to-backward ratio, 
F/B, is 1.17 +0.04, 1.06 +0.09, 1.08 +0.04, and 
1.13 +0.03 for distribution (c), (d), (e), and (f), re- 
spectively. 
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and 10%, respectively. Competing reactions, 
p+p>2n-+2n*+n7°, (3) 


occur with frequencies of about 32%, 7%, and1% 
for n=2, 3, and 4, respectively. 

The dynamical fit for each event was done by 
using an IBM-709 computer program (KICK).® 
The probability of a background reaction 7~+p—4 
prongs being wrongly identified as either Reaction 
(1) or (2) is negligibly small. When the prongs are 
wrongly transformed into the pp c.m. system as 
pions, there is much too little visible energy and, 
consequently, several missing 7°’s are required 
in order to conserve both momentum and energy. 
Secondly, the pion contamination in the p beam 
was 36%, and the pion cross section for produc- 
ing 4-prong stars small, namely, about 2.5 mb 
compared with 15 mb for )+p-— 4 prongs. 

Reaction D+ p> 27~+27*+7°. We identified 674 
cases of Reaction (1). These events satisfied the 
following conditions: (a) The x? for the hypothesis 
(1) was <9; (b) the value of the missing mass of 
4 prongs, interpreted as pions, was within 2.5 
standard deviations of 135 Mev; (c) the fitted p 
momentum was > 640 Mev/c; and (d) there was 
no fit for hypothesis (2). Of this sample, 410 
events (60%) had x? <2.5 and a missing mass < 250 
Mev. The still-unresolved systematic errors of 
our measuring system produce an extended tail 
in the x? distribution, so that the cutoff value of 
x? has to be increased generally by a factor of 
two or three. The angular and momentum distri- 
butions of the pions did not show any substantial 
difference between the samples with y?< 2.5 and 
x?<9. Our confidence is increased by the fact 
that the momentum distribution of 7s agrees 
very well with that of charged 7’s. 

The c.m. angular distribution for Reaction (1) 
is shown separately for 7~ and m* in Figs. 1(a) 
and 1(b). The 7~ distribution displays a forward 
peak, whereas that of the 7* is backward, contrary 
to the statistical model. As expected on the basis 
of C and CP,’ the distribution of 7* appears to be 
equal, within statistics, to the reflection of the 
angular distribution of the 7~’s. To analyze the 
shape of the distribution, we have reflected the 
m distribution about cos@=0 and added it to the 
7 distribution [Fig. 1(c)]. We have found it con- 
venient (without inquiring into its physical mean- 
ing) to express our distribution in terms of the 
expansion 2(a;P;(cos@), where P; are Legendre 
polynomials and the coefficients are (in units of 
numbers of pions): a,=135.3+2.2, a,=24.1+3.9, 
a,=18.2+5.1, a,=24.4+6.1, a,=19.4+7.1, a, 
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=10.0+ 8.9, a,=-14.3+8.9, and a,=-10.2+10,1. 
Note that the sharpness of the distribution makes 
the fourth order polynomial at least essential for 
describing the observed effect. The angular dis- 
tribution of the 7s does not show statistically 
significant deviations from isotropy. The momen- 
tum distribution of all charged pions agrees well 
with the phase space distribution, which predicts 

an average momentum (K),,=439 Mev/c. How- 
ever, for 7~ (1*) emitted in forward (backward) 
cone with |cos@| > 0.70, we observe (K),,=451+10 
Mev/c, whereas for the nonforward pions it is 

427+ 7 Mev/c. For those events that have simul- 
taneously a 7~ (7*) in the forward and a 7* (n~) in 
the backward cone, (K),y=495+15 Mev/c (482+ 26). 
It is interesting to note that the hydrodynamic mod- 
el of multiple-pion production in nucleon-nucleon 
collisions® predicts an increase of pion momenta 

at small angles. 

To determine whether the like pions (1+z+ or 
17”) have a tendency to be emitted in pairs in 
the forward peak [Fig. 1(c)], we have counted the 
number of Reaction (1) in which both 7~’s have 
cosé > 0.80 (< -0.80 for m*z*), If the probability 
of finding the second 7~ in this angular interval 
were independent of finding the first 7~ in this 
interval, then the probability that both 7™ mesons 
would be found there would be 1.9+0.1%.° We 
observe simultaneous forward (backward) emis- 
sion of pairs of 7~ (7*) in 15 cases of Reaction (1), 
or 2.2+ 0.6%. 

It has been reported that the angles between the 
like pions produced in annihilation are smaller on 
the average than the values predicted by the phase 
space.’° Our data confirm this observation. The 
phase space predicts a value of 1.79 for the ratio 
y= (number of pion pairs with angle > 90°):(num- 
ber of pion pairs with angle < 90°), whereas we 
observe y=1.56+ 0.08 for like pions. To find out 
whether or not this is caused, in part, by our 
angular distribution characterized by the sharp 
peak, we have evaluated y for only those pairs in 
which at least one 7~ has cosé@ = 0.80 (< -0.80 for 
tT), and obtained 7 =1.38+0.13. For comparison, 
we have also evaluated y for the pairs in which at 
least one 7~ has cos@ < -0,80 (> 0.80 for 7*) and 
obtained y = 2.05+ 0.25. These two values of y 
differ by 2.7 standard deviations and indeed indi- 
cate that the angular correlations between like 
pions are influenced by the sharp peak in the for- 
ward direction [Fig. 1(c)]. It will not be possible, 
however, to obtain a quantitative answer as to 
the importance of the Bose effect! and/or 1-7 
interaction without further extensive studies. 


Reaction )+p+27~+2n*. On the basis of x? 
and the missing mass, 136 4-prong events have 
been identified as Reaction (2), This is a 4-con- 
straint case, and the maximum acceptable ,? was 
taken to be 36; about 65% of these events have 
x?<10. The angular distribution of 7~ (1+) shows 
a forward (backward) peak [Fig. 1(d)]. 

Reaction p+p + 27~+2n*t+nn°, The remaining 
810 events, which did not fit either Reaction (1) 
or (2), were assumed to be Reaction (3). How- 
ever, this sample contains an estimated 16% of 
annihilations into K mesons, Reactions (4) to (7) 
(see below), in which K mesons have not decayed 
in the chamber. Because we transform our events 
into the c.m. system by ascribing pion mass to all 
charged prongs, the resultant c.m. angle will be 
wrong for these events. Also, this sample con- 
tains 4 prongs produced by 7~ in the antiproton 
beam, whose contribution could be as high as 5%. 
In spite of this background, the angular distri- 
bution of 7~ and m* exhibits features similar to 
the ones of Reactions (1) and (2) [Fig. 1(e)]. The 
average momentum is 369+ 9 and 332+ 5 Mev/c 
for forward and nonforward pions, respectively. 

(ii) Reactions j+p+K+K+nz. Of the 287 KKnz 
events, 135 were V-(2 prong) or V-(4 prong) events 
having a K*K® or K“K° pair from the reactions 











p+p+K* +17+K°+n,7°, (4) 

p+p+K* +1*+207+K°+n,7°, (5) 

p+p~+K~+1*+K°+n,7°, (6) 
and 

p+p +K~+07+ 20*+K°+n,r°, (7) 


As can be seen from the nature of Reactions 
(4) to (7), this sample contains a variety of types 
of events from the point of view of the fitting pro- 
cedure, with the number of constraints 1, 3, and 
4, Clearly each of these requires a separate set 
of criteria for identification. We shall give only 
a general outline: (a) The x? distribution for each 
type of event agreed well with the expected distri- 
bution when it was divided by 3; (b) the missing 
mass distributions were used to determine the 
number of neutral particles (™™s, or K° plus 7°’s) 
associated with each event; and (c) the data were 
examined and found to be resaonably free of bi- 
ases by looking at azimuthal distributions of events 
in the bubble chamber, their spatial distribution 
along the chamber, and the distribution of decay 


times for K° (or K°). The events were not gener- 
ally ambiguous as to the number of missing neu- 
tral particles. For most of the events, only one 
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of the K mesons was observed to decay within the 
chamber. Some events are ambiguous as to which 
particle is the second K meson. The detailed pro- 
cedure for the identification and disposition of the 
ambiguities of these events is given elsewhere. 
The angular distribution of both K~ and 7~ (K* and 
7) has exhibited a forward (backward) peak in a 
fashion similar to the ones of Reactions (1) to (3), 
as shown in the histograms in Figs. 2(a) and 2(b). 
Our angular distribution indicates a possibility 
that two different annihilation mechanisms are 
operative: the statistical one, which is dominant 
and produces an isotropic angular distribution; 
and a dynamical mechanism responsible for the 
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FIG, 2. Histograms of (a) combined angular distri- 
butions of K~ and reflected K* from Reactions (4) to 
(7); and (b) 7 and reflected 7* from Reactions (4) to 
(7). The F/B ratio is 1.37 +0.18 and 1.24 +0.15 for 
distributions (a) and (b), respectively. 
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forward peak. Whereas the observed forward- 
to-backward ratio F/B is consistent with the mod- 
el of Koba and Takeda,’ it is not clear whether 
the fact that kaons show similar angular distri- 
bution features is consistent with that model, be- 
cause it has never been extended in all its man- 
ifestations for the annihilations in flight. 

The very recent data on antineutron-proton an- 
nihilations at our momentum indicate a possible 
excess of 7* in the direction of the proton." 

It is a pleasure to acknowledge our particular 
indebtedness to L. W. Alvarez for his guidance 
and discussions, and to A. H. Rosenfeld for in- 
valuable aid in the concept of the data analysis 
for this work. Thanks are due to J. Button, 

P. Eberhard, A. Hussain, J. Lanutti, G. Lynch, 
N. Xuong, C. Tate, and C. Rindfleisch for their 
help in various stages of the experiment. 





*This work was done under the auspices of the U. S. 
Atomic Energy Commission. 

'E, Fermi, Progr. Theoret. Phys. (Kyoto) 5, 570 
(1950). For review on annihilation see E. Segré, Ann. 
Rev. Nuclear Sci. 8, 127 (1958). See also: P. Srivast- 
ava and G., Sudarshan, Phys. Rev. 110, 765 (1958); 

I, Pomeranchuk, Doklady Akad. Nauk S.S.S.R. 78, 
889 (1951); F. Cerulus, Nuovo cimento 14, 827 (1959); 
S. Frautschi, Progr. Theoret. Phys. (Kyoto) 22, 15 
(1959); B. Desai, Phys. Rev. 119, 1390 (1960); and 
G. Pinsky, G. Sudarshan, and K. Manhantappa, 
Proceedings of the 1960 International Conference on 
High-Energy Physics at Rochester (Interscience Pub- 
lishers, Inc., New York, 1960), p. 173. 

27. Koba and G, Takeda, Progr. Theoret. Phys. 
(Kyoto) 19, 269 (1958). This is mainly a model for 
annihilations at rest and does not discuss dynamic ef- 
fects. 

3G. Goldhaber, talk at CERN, November, 1960 (un- 
published); also, S. Goldhaber, Bull. Am. Phys. Sec. 
4, 417 (1959). See also reference 13. 

4G, Lynch, J. Button, P. Eberhard, G. Kalbfleisch, 
J. Lannutti, B. Maglié, M. Stevenson, and N. Xuong, 
Bull. Am. Phys. Soc. 6, 40 (1961); also, G. R. Lynch, 
Revs. Modern Phys. 33, 395 (1961). The background 
problem is now being investigated, because these events 
also fit reaction p+ p— 7* +27 +7°, 

5J, Button, P. Eberhard, G. Kalbfleisch, J. Lannutti, 











G. Lynch, B. Maglié, M. L. Stevenson, and N. H. 
Xuong, Phys. Rev. 121, 1788 (1961). Forward peaking 
is also observed in A emission from reaction p+p~A** 

®Reference Manual for KICK IBM Program, edited by 
A. H. Rosenfeld {Lawrence Radiation Laboratory Report 
UCRL-9099, 1961 (unpublished)]. See also A. H. Ros- 
enfeld and J. N. Snyder, Lawrence Radiation Labora- 
tory Report UCRL-9098, 1960 [Rev. Sci. Instr. (to be 
published)]. 

"A. Pais, Phys. Rev. Letters 3, 242 (1959). 


c OF © wt <a mic 


1 
1 
1 








rd- 
ie mod- 
ther 
stri- 
1, be- 
nan- 


yn an- 
sible 


ular 
nce 


‘sis 


mech, 
heir 


U.S. 


570 

, Ann, 
srivast- 
8); 

78, 
1959); 

, 15 
ind 


on 


Pub- 


S. 
or 
c ef- 


(un- 
Sec. 


isch, 
10ng, 
Lynch, 
ound 
events 


annutti, 
H. 


VoLUME 7, NUMBER 4 


PHYSICAL REVIEW LETTERS 





AucusT 15, 1961 





85, Z. Belenskij and L. D. Landau, Nuovo cimento 
3, 15 (1956). For the predictions of other models see: 
S. Takagi, Progr. Theoret. Phys. (Kyoto) 7, 123 (1952): 
W. Heisenberg, Z. Physik 133, 65 (1952); H. Bhabha, 
Proc. Roy. Soc. (London) A219, 293 (1953); and G. Coc- 
coni, Phys. Rev. 111, 1699 (1958). Charge distributions 
of the secondaries have not been considered in any of the 
ultrahigh-energy N-N interaction models. 

*This is strictly true only when the number of pions 
emitted in an annihilation becomes infinite. 

1G. Goldhaber, W. Fowler, S. Goldhaber, T. Hoang, 
T. Kalogeropoulos, and W. Powell, Phys. Rev. Letters 
3, 181 (1959). 


'G, Goldhaber, S. Goldhaber, W. Lee, and A. Pais, 
Phys. Rev. 120, 300 (1960). 

2G, R. Kalbfleisch, Ph. D. thesis, Lawrence Radia- 
tion Laboratory Report UCRL-9597, 1961 (unpublished). 

'3Hartmut Pilkuhn, University of Stockholm (private 
communication). 

4c. K, Hinrichs, Ph, D. thesis, Lawrence Radiation 
Laboratory Report UCRL- 9589, 1961 (unpublished). 
Also, it has been communicated to us recently that the 
data at 1.05 Bev/c in propane (reference 11) tend to 
agree with our effect, with the ratio F/B =1.18 +0.10 
[G. Goldhaber, University of California, Berkeley 
(private communication)]. 





SCATTERING OF Bev ELECTRONS BY HYDROGEN AND DEUTERIUM* 


R. M. Littauer, H. F. Schopper,f and R. R. Wilson 
Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 
(Received July 21, 1961) 


At Cornell University, using the 1.3-Bev elec- 
tron accelerator, measurements of the scatter- 
ing of electrons by hydrogen and deuterium have 
been extended to new energies and angles. The 
cross section at 90° has been measured using 
the same equipment and procedure that were 
previously employed’? at 66° and 112°. A new 
quadrupole magnet spectrometer of large solid 
angle has been placed at 45° or at 135° so that 
measurements could be made at two angles si- 
multaneously. In Table I are shown the new re- 
sults at 45°, 90°, and 135°, while all of the Cor- 
nell measurements are summarized in Fig. 1. 


The hydrogen cross section, do/d2, was de- 
termined in a straightforward manner. We have 
used the Schwinger correction for radiation. For 
the deuteron, the differential cross section d*oa/ 
dE'd& was measured at the energy where the 
electrons scattered from protons give an elastic 
peak. Then the total deuteron cross section, 
shown in Fig. 1, was calculated from this deu- 
teron peak cross section using the impulse ap- 
proximation as given by Goldberg.* Most recent- 
ly, Durand* derived Goldberg’s formula more 
rigorously and showed that it is accurate to with- 
in a few percent. 


Table I. Differential cross sections for elastic scattering of electrons of energy Ey by protons, and peak cross 


sections for scattering by deuterons. All parameters in the laboratory frame. Schwinger correction applied. 
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Ey Proton (do/dQ) x 10°? (em*/sr) Deuteron (d*0/dE’d) x10 (cm?/Mev sr) 
(Mev) 45° 90° 135° 45° 90° 135° 

317 6.04 +0,24 1.21 +0,08 2.14 +0,14 0.623 +0.080 
387 92.5 +4.0 3.61 +0.15 0.790 +0.031 31.2 +2.0 1.29 +0.05 0.307 +0.026 
407 81.0 +1.9 31.2 +0.7 

465 1.91 +0,06 0.296 +0.019 0.635 +0.04 0.130 +0.011 
552 37.5 +1.1 10.9 +0.3 

565 0.815 +0.08 ese 0.330 +0.03 eee 

600 0.99 +0.03 0.178 +0.008 0.333 +0,015 0.0638 +0,004 
664 22.0 +0.7 5.88 +0.3 

720 0.388 +0,017 0.0732 +0.005 0.132 +0.010 0.0248 +0.0032 
800 10.4 +0.5 2.68 +0.15 

836 0.212 +0.01 0.0302 +0.005 0.0448 +0.003 0.0121 +0.0026 
941 6.63 +0.3 1.39 +0.06 

974 0.080 +0,014 0.0272 +0.004 0.0147 + 0.003 0.0083 +0.0024 
1050 4.68 +0.2 0.76 +0.04 

1136 0.0362 +0.0028 0.0117 +0.004 0.0107 +0.003 0.0033 +0.0015 
1166 3,24 +0.17 0.61 20.06 
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Differential elastic scattering cross section for the proton (a) and the deuteron (b), as a function of 


incident electron energy and laboratory scattering angle. Full curves computed with core model, dashed curves 
according to Bergia et al. See R. M. Littauer et al., following Letter [Phys. Rev. Letters 7, 144 (1961)]. 


Figure 2 illustrates the determination of the electric form factors F, and the magnetic form factors 
F, by the intersecting ellipse method.* For a given value of q*, the square of the momentum-energy 
transfer, one ellipse in the F,-F, plane is determined by a measurement of the cross section at a par- 
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Proton 
' q*=7,41-2 





Neutron 
' q*= 103? 


FIG. 2. Determination of F; 
and F, by the method of inter- 
secting ellipses. Examples for 
the proton at q?=7.4 f? (a), 15.0 
f~? (b), and 24.7 f~* (c); for the 
neutron at g?=10.3 f-* (d). In 
each case the abscissa is F,, the 
ordinate F,/p. 
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FIG, 3, Electric and magnetic form factors for the 
proton and neutron. The points are obtained by the in- 
tersecting-ellipse method; thus, errors for F; and F, 
are not independent. Typical acceptable deviations are 
indicated, in correlated senses, by the solid and broken 
arrows. The solid curve shows a fit to the data by the 
core model, the dashed curve by the Clementel-Villi 
form, 


ticular angle. Other measurements at different 
angles but at the same gq’ determine other ellip- 
ses. (The form factors are presumed to be func- 
tions of qg? only.) The intersections® give specific 
values of F’, and F,, which have been plotted 
against g? in Fig. 3. 

In our previous measurements of the scatter - 
ing from deuterium we had results at but one 
angle; hence, in analyzing our data we neces- 
sarily made the assumption that F 9, = -F 9». 
This assumption, of course, is now no longer 
made in finding the results shown in Fig. 3; we 





find, in agreement with the recent Stanford meas- 
urements,’ that -F'9, becomes larger than F»9 
especially at high values of g?, where F2p is 
small or might even go negative while -Fo, is 
still positive. These new values of F9, do not 
seriously change our previously deduced values 
of Fip (by less than the statistical error), but 
they do indicate that it is necessary to assign a 
small negative magnetic moment to the core of 
the proton and neutron which gives rise toa 
small isoscalar part of the magnetic form fac- 
tors. 

Our values of F4y,/uy are positive and between 
0.1 and 0.2, in agreement with the earlier Stan- 
ford results’ calculated with a modified Jankus 
theory.®»> This disagrees, however, with the 
reinterpretation of the Stanford results by Du- 
rand,* which yields values of Fy,/uy, that are 
small or may even go negative. 





*Supported in part by joint contract of Office of Naval 
Research and the U. S. Atomic Energy Commission. 

TOn leave from the Technische Hochschule Karlsruhe, 
Germany. 
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STRUCTURE OF THE PROTON AND NEUTRON* 


R. M. Littauer, H. F. Schopper,f and R. R. Wilson 
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We have previously analyzed our data’ to show 
that the proton and neutron have a common core 
of positive charge of about 0.35e, whose radius 
might be about 0.2 fermi. This core is sur- 
rounded by a meson cloud of radius 0.8 f, with 
an average charge which changes from +0.5e 
for the proton to -0.5e for the neutron. Addi- 
tionally, we found it necessary to postulate 
another cloud, the same for the neutron and 
proton, with a positive charge of 0.15e and an 
rms radius of about 1.4 f, thus extending out 
beyond the meson cloud. These clouds of charge 
take on meaning in the simple A model of the 
nucleon put forward by Bergia, Stanghellini, 
Fubini, and Villi,? referred to herein as BSFV. 
Their model is based on dispersion theory and 
on a two-pion interaction similar to that used 
by Frazer and Fulco.’ In this model they identify 
the meson cloud with a T=1, J=1 two-pion reso- 
nant state while the extended cloud might corre- 
spond to a three-meson resonant state, most 
probably with T=0 and J=1. Our new measure- 
ments,* when examined in the light of the above 
model of the nucleon, imply that the resonant 
energy of the two-pion state is 4.0m, and that 
of the three-pion state is 2.9m,. 

Let us now analyze our form factors into iso- 
scalar and isovector partial form factors’: 


=F ° ¥ = - 
Fg ig aye on 18? we 


Fay*FagtF ays Fay *Fas-F ay: (1) 


In making our simple core model analysis, we 
got a vanishing rms radius for the neutron charge 
distribution by resolving the isoscalar form fac- 
tor Fj, of Eq. (1) further into two terms, i.e., 
F45=F 15°’ +F 15°, where we associated 
F15°°F®, the term with the smaller radius, with 
the core of the nucleon, and Fis“, with an ex- 
tended cloud. Our new results can best be fitted 
by partial form factors corresponding to exponen- 
tial density distributions. Partial form factors 
of Gaussian or other shape without a singularity 
give a worse fit. We write therefore 
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22 -2 
1S os s,core q /12) 


cl 2 2 -2 
eae” (ite, q /12) , 


= 2 -2 (2) 
Foy ef +? q°/12)~, 
where the a’s are the rms radii of the respective 
charge distributions. In order to obtain the right 
total charges and rms radii for the proton and 
neutron, the following conditions must hold: 


core cl 
e +e 
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e +e -e =0, 
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core 2 cl 2 2 2 
e a + a - A 
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core 2 cl 2 2 2 
e a +@ a 
Ss s, core s Se vv n 


s 2 2 2 
Assuming that as core <« as cl and because a, 
=0, one obtains 


2 cl 2 2 (4) 


These conditions restrict the six adjustable pa- 
rameters so that only two are left to fit the ex- 
perimental data at higher q values, e.g., e,©F® 
and 4s, core- 

To describe our newly determined magnetic 
form factors, we now assign partial magnetic 
moments p,©°re, y Cl, yu, to the same core 
and clouds that we have used to describe the 
charge distributions and we write 

core 
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where bs core’ bs, c])» and b, are the rms radii 
of the partial form factors for the magnetic mo- 
ment and all magnetic moments are measured in 
nuclear magnetons (nm). 

In exactly the same manner as for the partial 
electric charges, we can infer the conditions 


core cl 
+p 


‘ 3 ia ih 1.793 nm, 


core cl 
+u -p =p =-1.913 nm, 
Ss Ss v n 


core 2 cl, 2 ~%" » 2 
s s,core *¥s "scl *Hy’y “Hp? 


core e cl, 2 b = b? 
s s,core ~"s "scl “Py ’y “Han * 


(6) 


From these one readily obtains with b, core’ 
«bs el”: 


yw =1, 
v 
core cl 
= = -0, 7 
M +H, (utu,)/2 0.060 nm, (7) 
cl 2 2 2 
au, .. cl =H 5? +u,>, . 
2u beau bub. (8) 


bp 


A consequence of these conditions is again that 
only two parameters, e.g., u,°°F© and by core? 


are left to fit the data at high q values, if by" and 


by? are obtained from the slope of the form fac- 
tors at small q values. 


Table I. Best-fit parameters for core model with ex- 


ponential density distributions. See text for definition 
of symbols. 








« =". os a =0.2 f 
Sal s,core 

e = 0.25 a =1,.13 f 
Ss s,cl 

e = 0.5 a =0.80f 
v v 

ae, -0,22 b undetermined 
Sal s,core 
= 0.16 b =1,30f 

s s,cl 

bu = 1,853 b =0.,89f 
v v 

a, = 0.80 f b =0.98 f 
a = 0 b =0.79f 








At large values of q’, the effect of the extensive 
isoscalar cloud will be negligible and we can fit 
the core parameters. A point moment of -0.19 
nm would give a satisfactory fit, but let us as- 
cribe to the core, arbitrarily, the same rms 
radius of 0.2 f which fits the charge core. In 
this case, the core moment must be -0.22 nm. 

In order to make the static moments correct, 
we must then ascribe a magnetic moment of 
+0.16 nm to the extended isoscalar cloud so that 
the total isoscalar moment = a(Hp +H) = -0.060 
nm, as required by (7). 

Table I summarizes the parameters we have 
found for the partial form factors; these are 
plotted in Fig. 1, and the corresponding spatial 
distributions are shown in Fig. 2. 

The qualitative result then of these new meas- 
urements is that we find it necessary to attribute 
a magnetic moment of negative sign to the charge 
core that our previous measurement had revealed. 
It is not unreasonable that the angular momentum 
of the meson cloud be rar es then necessarily the 
core will have a spin of -4, which might rather 























F; 
1O 
Br 
Br 
Fe Foy 
4+ 
Ss 
ae 
a Bs fcore 
° 4 2s aes 
- 1 I 1 
= 0 20 30 
FIG, 1. Partial form factors for F; and F,. Each 


form factor is resolved into isoscalar and isovector 
parts, whose sum and difference give, respectively, 
the neutron and proton form factors. The solid curves 
indicate the fit according to the core model, where the 
scalar partial form factors have been further split into 
terms corresponding to a core and an extended cloud, 
each of exponential distribution. The dashed curves in- 
dicate the best fit obtained by the Clementel-Villi form. 
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naturally give rise to the negative magnetic mo- 
ment—for example, by dissolution of the core 

into a K*A system in direct analogy to pion emis- 
sion. It is also perhaps significant that the radius 
of the extended isoscalar cloud of magnetic mo- 
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(l-a )+ 4 ’ 
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where one minus the constants a or 8 give the 
fraction of the charge or magnetic moment in a 
point core, and a, or a, give the corresponding 
fraction in a Yukawa meson cloud of mean square 
radius 6/q,* or 6/q,”. It appears that BSFV have, 
not unreasonably, assumed that the radii of the 
Yukawa clouds that appear in F, and F, are the 
same. They show that the radius of the Yukawa 


form factor is due 


to a T=1, J=1 two-pion resonant state with the 
resonance at q,, while the isoscalar radius will 
correspond to a three-meson resonant state, 
most probably with T=0, J=1, at q,. The six 
parameters in (9) and (10) are reduced to five 
because the experimental requirement that the 
neutron charge radius is zero leads to the rela- 


Table II. Parameters for best fit according to BSFV* 




















form. 
Cornell Stanford 
a 1.10 1.20 
v 
B,, 1.14 1.20 
a. 0.58 0.56 
B. -1,5 -3.0 
a, 0.85 f 0.77 f 
a. 1.16 f 1.13 f 
a, 0.88f 0.85 f 
b, 0.95 f 0.94f 
b 0.87f 0.76 f 
q,° 8.3 f= 16(m_/c)* 10 f-?= 19.6(m_/c)? 
q,” 4.4 f*=8.5(m_/c)? 4.7 f= 9(m_/c)? 








a 
See reference 2. 
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ment comes out so close to that of the isoscalar 
charge cloud—also true of the isovector clouds. 

BSFV have proposed aesthetically simple form 
factors for the nucleon based on the model already 
referred to above. In our notation their form fac- 


tors are 
v e 
F, -1.889{( -B) ‘aah (9) 
F,*=-0.060|(1-8 )+ °s 10 
2 = "Ve s 1+q°/q" , ( ) 





a /q*=a /q?*. (11) 
s Ss vv 


Table II gives values of the five independent pa- 
rameters that best fit our data plus the rms radii 
of the neutron and proton. The same parameters 
have also been determined at Stanford® and their 
values in our nomenclature are shown for com- 
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FIG. 2. Spatial distribution of charge and anomalous 
magnetic moment for proton and neutron, according to 
the core model. 
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parison; the general agreement is good. Slight 
discrepancies arise mainly because our experi- 
mental values for Fo, at high g* values are close 
to zero whereas an extrapolation of the Stanford 
form factors gives negative values. Our experi- 
ments indicate also a somewhat smaller differ- 
ence between Fo, and Foy. 

The experimental data are fit equally well by 
the BSFV form factors or by those following from 
our simple core model (see Fig. 1 and the curves 
in the preceding Letter‘). 

According to the interpretation of BSFV, our 
values of q, and qs would imply that the resonant 
energy of the two-pion state is 4.0m, and that of 
the three-pion state is 2.9m ,. 





*Supported in part by joint contract of Office of Naval 
Research and the U. S. Atomic Energy Commission. 

TOn leave from the Technische Hochschule Karlsruhe, 
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X-RAY YIELDS IN THE K AND L SERIES OF ». -MESONIC ATOMS* 
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(Received July 14, 1961) 


Since the experimental work of Stearns and 
Stearns,’ much theoretical speculation has been 
presented?~’ to account for the low yields of x- 
rays in the K and L series of light u.-mesonic 
atoms reported by these authors. Their results, 
reproduced in Figs. 1(a) and 1(b) (filled circles), 
require that the Auger rates® competing with radi- 
ative transitions be, respectively, ~300 times 
(K series) or ~30 times (L series) as fast as pre- 
dicted by conventional theory, or else that one 
postulate some hithertofore unknown competitive 
transition mechanism. The situation for the yields 
of 7-mesonic x rays of comparable energies is 
entirely similar,® but need not be discussed here. 

In view of this situation, we decided to remeas- 
ure the yields in the 7-mesonic K and L series, 
in particular for those light elements for which 
low (or even experimentally unobservable) yields 
were reported in reference 1. In view of our find- 
ings, we think it useful to present here our pre- 
liminary results. They are displayed in Figs. 1(a) 
and 1(b) (open circles). 

The experimental arrangement used by us is 
Shown in Fig. 2 and is essentially the same as 


that used by Stearns and Stearns.’ Some relevant 
differences are the following: (a) use of thinner 
(1/16 in.) plastic scintillators, wrapped in 1-mil 
Al foil, just before and just behind the x ray tar- 
get T, to improve the transmission of soft x rays 
and to decrease the background of carbon mesonic 
x rays; (b) use of an improved window (20 mils 
of Be) on the 1/16-in., 2-in.-diameter Nal de- 
tector employed for x rays up to 75 kev; (c) use 
of a longer effective resolving time (100 nsec vs 
the 50 nsec of Stearns and Stearns) for NalI-in- 
duced coincidences; (d) use of a 100-channel 
pulse-height analyzer (RIDL 34). 

The running conditions differed more markedly 
from those available to earlier workers in this 
field, in the following respects: (a) use of a u~ 
beam of low e~ content and of a Cerenkov counter 
(No. 3 in Fig. 2) for further electron rejection; 
(b) use of a good duty-cycle (~5) meson beam, 
obtained by means of a vibrating cyclotron target.’® 
These conditions enabled us to use rather slow 
Nal coincidences without generating large acci- 
dental backgrounds. As is well known, hard clip- 
ping of NaI pulses of small amplitude (from, say, 
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FIG. 1. (a) Observed yp K-series yields: open circles, 


this experiment; filled circles, reference 1. Curves 
are fits to nx =const Z‘/(Z*+C,). For this experiment, 
const=1.00, C,=45, while const=0.86, C,=450 for 


reference 1. 


With the same assumptions as used in 


reference 1, present data require an Auger transition 
probability about 30 x theoretical. Note that the sole 
evidence for fall of K yield with low Z rests on the Li 


K point. 


(b) Observed » L-series yields: open circles, 


this experiment; filled circles, reference 1. Curves 
are fits to my =const Zzi/(z4+ C,). For this experiment, 
const= 0.80, C,=3 10°, while const=0.97, C,=2 10! 


for reference 1. 


With the same assumptions as used 


in reference 1, present data require an Auger transi- 
tion probability about 4 x theoretical. 


25-kev photons) can easily lead to large time jitters and hence appreciable coincidence losses. The 
limiter (see Fig. 2) used by us was well cut off by 12-kev x rays from a Hg”™ source under the oper- 
ating conditions of the NaI detector; further, during actual runs a good plateau (<2% slope/100 v at 


2100 v photomultiplier voltage) was obtained with Li K x rays (19 kev). 


150 Mev/c 
(we) beam ———> 
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Vibrating Target—> 


2 3 45 6x 






















































FIG, 2. Experimental arrange- 
ment: 1, 2, 4, 5, and 6 are plas- 
tic scintillators, viz., 1=6*6 
x 3/8 in., 2=3*3%1/4 in., 4=5 
x5X1/8in., 5=3*3%1/16 in, 
6=4.5*4.5 x 1/16 in., 4=4%4 
x1in. H,O—Cerenkov counter, 
X = 2-in. -diam <1/16-in. Nal 
crystal, connected by logarith- 
mic light pipe to a 5-in. -diameter 
photomultiplier (10-stage EMI). 
A;, Ag, Aj absorbers; A, =2 in. 
graphite, A,=0.75 in. Cu, A; 
= 0.008 in. brass (shield for x 
rays from counter 5); 7, target 
(in Table I; L = limiter; CF =cath- 
ode follower; HP = distributed 
amplifier, gain 10; D = fast dis- 
criminator; G=5-msec gates; 
A=linear amplifier; PHA=100- 
channel analyzer. (1,2), (14536), 
and (py) are coincidence circuits 
with the indicated resolutions. 
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Table I. Experimental data used for yield evaluation. 








Element Target Series, Stopped” Calibration Counting” Detected“ _Kq_ Series yield © 
thickness Kg (Lg) Muons x 10° source efficiency X rays all K Muon 
(g/ cm?) (kev) (corrected) correction (%) 
Li Metal, 1.2 K, 18.7 35.1 18 kev, Am™! —s_ 9,091 1.700 0.70 0.53+0.05 
Be Metal, 1.5 K, 33.3 76.3 32 kev, Cs'9? 0.078 6.061 0.82 1.01+0.10 
Cc Graphite, 1.2 K, 75.0 18.9 73 kev, Hg? 0. 086 1.592 0.71 0.97+0.10 
Cc Graphite, 0.48 L, 13.9 44.6 12 kev, Hg? 0.067 0.475 0.63 0.16+0.05 
B Powder, 1.9% K, 52.1 52.0 52 kev, Tm'™ 0,093 5.475 0.77 1.12+0.17 
N 95 % (NHq)» K, 102 29.9 60 kev, Am”! 0.16 4.836 0.69 1.01+0.10 
+5% HO, 1.0° L, 18.9 34.6 18 kev, Am™! —s 0,077 0.889 0.63 0.34+0.03 
oO Ice, 1.0 L, 25 32.7 32 kev, Cs‘? 0.088 0.964 0.68  0.32+0.03 
F LiF, 1.0° i. = 23.2 32 kev, Cs‘? 0. 087 1.334 0.81 0.66 +0.10 
Mg Metal, 2.0 L, 56 36.6 60 kev, Am™! 9), 091 2.892 0.75 0.87+0.04 





“Contained in 3-inch diameter x 1/16-inch wall bakelite cylinder with 0.02-inch Lucite end windows. 
oF faction of muon stops in other than target material was measured. 
Includes empirically determined target, ac counter, and Nal casing attenuation, in- and out-scattering, and 


Nal efficiency (escape correction theoretical for Be K). 


a NT runs with 2-inch diameter x 1/16-inch Nal detector except N K with 3-inch diameter x 1/4-inch Nal detector; 
C K background was measured independently, C L background from counter No. 5 was filtered out with 0.008 inch 


brass. 
“indicated errors are increased over statistical to include 
in absolute calibration with sources. 


Table I summarizes the experimental data lead- 
ing to the yields plotted in Figs. 1(a) and 1(b). The 
targets were in general 3-in. diameter cylinders, 
matched in surface density for roughly equal at- 
tenuation (~10%) of the mesonic x rays under 
study, rather than for equal stopping power,’ 
and subtending approximately equal solid angles 
at the detector; the number of stopped particles 
was determined experimentally. The targets 
were calibrated in situ with uniform extended 
radioactive sources of calibrated intensity" yield- 
ing photons of energies close to those of the mes- 
onic x rays of interest. Further calibrations, 
mainly for in- and out-scattering, were performed 
in the laboratory by displacing extended sources 
through appropriately subdivided targets. The 
important escape correction (~35%) for the Be K 
radiation (33 kev) had to be taken from theory’? 
because no radioactive sources yielding photons 
as close to the iodine K edge (33.2 kev) were avail- 
able. Dead-time corrections to the data stored 
in the pulse-height analyzer were negligible (10 
analyses/sec, 100-ysec analysis time). 

As is evident from Figs. 1(a) and 1(b), our re- 
Sults do not appear to bear out those of reference 
1, Assuming that the mesons have a statistical 
distribution of 2/+1 in the n=14 level, their re- 
sults require an Auger rate of 5.85 x10'* sec™ 
(~300x “theoretical”) for the K series and an 


error in analyzing spectra for background and error 


Auger rate of 2.68x10"* sec (~30x “theoretical”) 
for the L series. In both cases the remaining 
“discrepancy” is an order of magnitude less seri- 
ous than suggested by reference 1. By making 
different, but not unreasonable assumptions about 
the initial populations,’* one can hope to reduce 
the extant “discrepancy” even further. 

It has occasionally been suggested that the re- 
sults of Stearns and Stearns’ could be attributed 
to a thin “dead layer” on the front face of their 
Nal(T1) crystal, or of such crystals in general. 
This hypothesis is not only at variance with gen- 
eral experience with NalI(T1l) detectors, but also 
fails to provide a satisfactory explanation. The 
effect of such a “dead layer” should be quite dif- 
ferent for x ray energies just below, and just 
above, the iodine K edge, i.e., for Be Ky (33.3 
kev) and F Ly (32 kev). The yields of reference 
1 show no great variation between these two x 
rays. 

Another suggestion frequently advanced is that 
the muon may be “trapped” in some high-n orbit. 
We briefly investigated this possibility by delay- 
ing the py input to the p.-y coincidence circut (see 
Fig. 2) by 100 nsec and looking for Li K x rays, 
the ones for which we find the lowest K yield. Not 
more than (1.0+ 0.5)% “delayed” x rays were 
found in this way. This observation is most nat- 
urally interpreted in terms of some residual jit- 
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ter of the NaI pulse, but can clearly not preclude 
the existence of a long-lived delayed component. 
In the case of the 7~ mesons, recent work at 
CERN" suggests, however, negligible trapping 
in Be. 

A repetition of the present .-mesonic yield 
measurements with a proportional counter is in 
progress, and an extension to pions is contem- 
plated. 
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